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§1. Introduction 


In another paper [1] we study a Riemannian space”) of dimension 2n 
admitting a Killing vector 4 such that 


(1.1) Vau' (Fa ur)=Gpa 


and obtain following results. 
This Riemannian space V2” becomes an almost Kahlerian space by virtue 
of the tensor?) 


(1.2) F34=V3u4 . 


At any point where the Killing vector does not vanish we can take 


a coordinate system (&4) such that 


(1.3) u4=6 , Gsa=0 
and 
(1.4) Goo=&" , Gos=Gie=0 . 


Such a coordinate system is called a favourable coordinate system. For 


this coordinate system Gaza satisfy the equations 


i} 
‘e000 — = 79 pose - ; 
(1.5) GonZo, Greate 
(1.6) E°O.Gor=Gor , 


(1.7) GG Gop} 5 (Gu Go G4 =Gi EY "GG | 


We define fi and gji by 


1) The fundamental tensor of this space V2" is denoted by Gra and the coordinates 


by £9, £1,-.-, g2-2, ge, We use indices as follows, 
AUB Cl, 8) Ty =0, 1, 2a 2, 
EC 
DE PY =0, 1,・…・, 2n-2, 
h, 4,9, k,l, m = en 


2) Only local properties are considered. 
3) We use the symbol F's4 or fj% to denote an almost complex structure. 
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Goi 


7 — (fo)\-1 i 
(1.8) Kk Goo (E°)-!Go 
and 

Gii GojGoi 
Gy) Set Guan 


Then fi do not contain the variable &”, while gj;i depend in general on &°. 
We obtain a family {Vm-2(E) Es D} of Riemannian spaces of dimension 2”—2, 
with the fundamental tensor g;i and such that each point is denoted by the 
coordinates &!,---, &2"-1, while to each value of &°€ D corresponds each member 
of the family.” V2"-2(E0) becomes an almost Kahlerian space by virtue of 
the tensor 


(1.10) f#=fag" , i= fi 


Conversely, we consider a family of Riemannian spaces Me”) (&°€ D) of 
dimension 27—2 such that each point in a member Me) of the family is 
denoted by the coordinates &',---,€"-2 and such that the fundamental tensor 
is given by the components gji(&!,---, €"-?; €°). Moreover we assume that 
each Me) admits a vector field such that the components f; do not involve 
the variable €° when the coordinate system (é',---,&"-2) is chosen suitably 
and such that Me) becomes an almost Kahlerian space by virtue of the 
structure fj'=fixg™, fii=3(0; fi—0:f;). Then a Riemannian space of dimension 
2n where a point is denoted by the coordinates &°, &',---, €2"-2, © and such that 
the fundamental tensor Gza is given by 


(1.11) { Go=& ’ Goo=Gin=0 ? Gr = (4é~)-1 : 


Gor=Er"fe,  Gr=eECOR+SA) 
admits a Killing vector with components w4=d 94 and this vector wu satisfies 
(1.1). Hence this Riemannian space V*" is an almost Kahlerian space with 
Fi4=Vaud and its V"-(€) is M(E). 

We also have the following theorems. 

A necessary and sufficient condition for an almost Kahlerian space V2" 
with Fsi=Vasu‘® to be a Kahlerian space is that Vn-? be Kahlerian spaces 
and that the hypersurfaces uur=const be totally umbilical. 

A necessary and sufficient condition that, in an almost Kahlerian space 
V*" with F34=Ve2u4, hypersurfaces u’ur=const be totally umbilical is that 


(1.12) £°0.G = Gii 


4) D is a domain of positive numbers. 


5) It is evident that, if Gaa4 and F'34=p, u4 determine an almost Kahlerian structure, 
then u4 is a Killing vector. 


a 
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hold for any favourable coordinate system. 


§2. Special Kihlerian spaces® 


Since (1.12) is equivalent to that the functions gji do not contain €°, we 
obtain the 

THEOREM 2.1. A necessary and sufficient condition that an almost Kahlerian 
Space Vr with the structure Gra and F34=Vzu4 be a Kahlerian space is that 
the spaces V*"-*(E~) are isomorphic to a Kahlerian space Vo2"-?.” 

We calculate the components F34 in a favourable coordinate system. We 
get 


0 =G"l00, 0]-+G00, N=0, 


{ Gli0, 0 + 6170, =F GO;Gu—2Gu)=—S'fn=—Sithe 
US Pe 00 ol BAS 00 1 or Gor _ 1 
{of =6 [oo0, 0] + G"[o00, = + Gr =o, 

? \ Go0, 0] +G*{00, N=0 
{ o0 | =@i"t00, 01+ (00, 7]=0 , 
{ : been 70, 0) + G0, N=+G"@sGu—-IGos)=9''fin= rae 
t | pio il seat iy ported te i Goi _¢ 
{cof =2 [oo0, 0] + Gloo0, 1] =3-G* + = Gt" =0, 

co eee ee Soe 
{oc [00, oo]=—3G 26° , 

fore) 3 G-? Gos 

YES pe eee 

1; } [70, 00] ake fi 


fh cca co]=0 


by using Gza=0) (1.11) or (1.4), (1.5), (1.6), (1.8), (1.9), and (1.10). As we 


; A A 
have Fit=dout+| he =f Bob we get 


Fj;°=0, F;'=0, Fjr=—2e , 
(2.1) | OP aah Fe ee Bia ij ts, aah = 
1 ie : 
oe ’ Fs=0 » Fes0 . 
F 5 


6) In the present paper also we consider local properties. 
7) An almost Kahlerian space V, is isomorphic to an almost Kahlerian space V2 
when an isometric transformation V,—>V2 conserves the almost complex structure. 
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Let us consider a Riemannian space V of dimension 2n—2 and denote its 
fundamental tensor by gji. We assume that it admits a Kahlerian structure 
fj* such that fji=3(0;f.—d:f;)g" for some vector field fj. Then we have 


* = 
fifit=—0,  Vefii=0 
. 
where Vi; denotes the covariant differentiation with respect to the Christoffel 
* 
symbols 20 formed from g;i. From this space V we construct a space 
jt 


{Vx&xé|E°E R, €°>0} and try to introduce an almost Kahlerian structure 
Gea, F's“ by (1.11) and (2.1). 

As we have F3%Fs4=—6z4 from (2.1), Fs4 is an almost complex structure. 
As we have Gsz4=0, u4=d04 is a Killing vector. Moreover, we have determined 
F34 in such a way that we may have Fz4=Vzu4. Hence we have an almost 
Kahlerian structure. This becomes a Kahlerian structure according to Theorem 
2.1. Thus we obtain the 

THEOREM 2.2. Let V be a space of dimension 2n—2 admitting a Kahlerian 
structure gyi, fj*=3(0;fr—Onf;)g"*. Then the space {VX&xé&|—&ER, E°>0} 
admits a Kahlerian structure Gea, Fz4, where Gea is given by (1.11) and F's“ 
by (2.1). 

From this theorem we get a way to construct a Kahlerian space of 
dimension 2” from a Kahlerian space of dimension 2n—2. Thus we can obtain 
a sequence of special Kahlerian spaces V2, V*, V°,--- starting from a Kahlerian 
space V2?. 


§3. A special Kahlerian space V? 


Since V? is the first space in the sequence, we can take various spaces as 
V*. But if we consider that V? is also a Kahlerian space where F34=Veu4 
(A, B=0, oo), then this space is obtained formally from an imaginary zero- 
dimensional space. Thus V* is a space where a point is denoted by coordinates 
€°, & and the structure is given by 


(3.1) Cases CaeCn=6. Guo 
(3.2) Fe°=0, Fyr=—2&, Fr= F:°=0. 
We get 


Ko0=G{[00, 0][eccec, 0]—[000, O][0co, 0]} 
+G"{[00, co][coco, co]—[0co, co][Qoo, oo]} 
1 1 1 


— ———_(500 + —G»r~ 


4 4° 4 


0, 
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hence V? is flat. If we call a Kahlerian space such that F34=Vnu4 a trivial 
space when its curvature tensor vanishes, this V2 is a trivial space. 

If we use a rectangular coordinate system (x, y) and denote a Killing 
vector by (uw, v), we get from 0.u=0, 0,u+0.v=0, dyv=0 


u=—cyt+a, v=cx“xt+b. 


If c=0, the magnitude of the Killing vector would be constant. Hence we 
must think c#0. Choosing the origin suitably we can write w=—cy, v=cx. 
As the value of c does not concern us, we put c=1 and obtain 

u=—Y, =i 


The relation between (&°, &°) and (a, y) is obtained as follows. As (é)!/2 
is the magnitude of the Killing vector, we have °=x2?+y’=r?. We substitute 


this into 
B A Qa B 0 A 0E8 0 A 
Cus o=1 : Gus =0 < GG oo) 

and get 

0& +y UE 

Be ey? Cy OY ” 
hence 

S&=+0+k 


where @ is the argument of the polar coordinates (7, 0). The double sign is 


removed on account of w4=06o4 or 


0&° 0f° 
——b + ——U = 1 
Ox e oy 
and we get 
"=0 


for we can put k=0. 
Consequently we have 


=F, sr’ 


as the relation between the favourable coordinates and the polar coordinates. 
Though the point s=0, y=0 in rectangular coordinates is a singular point 
of the favourable coordinates, it is not a singular point of the tensor field 


F34, for its components are given by 


(Fad=( 7 


in rectangular coordinates. 
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Thus V? is a flat space locally, but we have various cases concerning its 
global form. For example we can consider that V? is a 2-dimensional euclidean 
space R2 We can also consider the case where V® is a space obtained by 
removing a point from Rs or one of its covering spaces if we do not assume 
completeness. According to such circumstances 9 must be restricted by 029 
<2z or 0<0<2mz. Another example of V? which is not complete is obtained 
by restricting r=(é")'/2 and 6=£ by 0<a<r<b and 00S and identifying 
the point (r, w) with the point (7, 0). 


§4. A special Kahlerian space V‘ 


Let us construct a V* from a V? such as is obtained in §3. For this 
purpose we must change the old indices 0, co in V? for new indices 1, 2 in 
V‘ and introduce indices 0, co anew in V‘. Gaza, u4 of V? changes into gji, 
f* of V‘. Then we have f!=1, f2=0, gu=&, gis=g2=0, g22—(46*)*, and 
consequently, fi=&, f=0. Hence the fundamental tensor Gza of V* is as 


follows, 
fo plead 0 0 
2 £co 2\ £2 Eco 
(chy で 0 0 
0 0 (4E i ed 0 
0 0 0 (4€=)-1 


Since we have from (3.2) 
ft=0.,4 fui=—2e.” Ft=2y1, fst=0, 


F34 has the form 


0 0 0 _2¢= 
0 0 282  _pg2g~ 
F34)= 
oe —} (282) 0 0 
(26)-1 0 0 0 


Let us change the real coordinates for the complex coordinates (x, Y, £, ¥). 
Changing the indices 0, 1, 2, oo for 1, 2, 3, 4 and changing the letter G 
for g, we obtain 


ct ee 0 0 
Gia | aati dees ee 0 
: 0 0 (4€8)-184 0 


0 0 0 (4€4)-1 
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esta?) —(E4)=1 0 0 
+ —(&4)-1 (Ge 0 0 
(GE 
0 0 (&4)-14€8 0 
0 0 0 4é4 : 

0 0 0 —2¢4 

(Fj*)= 0 0 — 23° soe: 
。 —4 (2€8)-1 0 0 
(264)-1 0 0 0 


Complex coordinates are obtained by solving the differential equations [2] 


(4.1) Bede 20), 
(4.2) Crd 
where 


Bi GAS CH= GA+IF . 
As (4.1) is equivalent to 


Vee 1 
RE AE 
de de = E05; 


; 1 
dé?+—d&=0, 
1067+ DE で 

we obtain for example 


z= ーー log &, 


ya + log &3 . 


Hence we have 


1_t—@ = UseY 
E 24 ? € 21 > 


=exp (y+Y) » 
4=exp {exp (y ギ ひ ) 十 〆 十 %} . 


Calculating the components of the fundamental tensor in this complex 


coordinate system, we obtain 
ds?=exp {exp (y+ y)+u+x}dudx 
+5 exp {exp (y+y)t+ta+z+y+y)}(dxdy+ dudy) 


+{1+exp (y+y)} exp {exp (yt+y)+a+Et+y +y}dydy . 
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In a Kahlerian space there is locally a real analytic function U(x*, x”) such 
that the components gaz of the fundamental tensor in complex coordinates 


(gl, 0", ates, a”) satisfy the equation 
ehh 
0 = RIE 


In our V: we have 
Yate exp {exp (y+y)+a2+Z} . 


In [1] we obtain the following formula showing a relation between the 
curvature tensor Kyors of V2" and the curvature tensor ee of ye 3 
(6°)-*Knjsa—_Kensen=—guagse+ gingua—( fen fs— fn fui) +2 fs fin 
+f0 ee IT. iFint {Ws fus— iV nt « 
In the present case yin is a Kahlerian space, so that we obtain 
(4.3) Kua {Kein = OTA dated : 


We can use this formula in calculating the curvature tensor of V‘ and obtain 
for example Kss23 =l. Hence V‘ is not a trivial space. 
(4.3) also proves that the special Kahlerian spaces V‘, V°,--- are not 


trivial spaces. 
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(Received June 30, 1961) 


Synopsis 


Ferromagnetic resonance experiment on single crystals of NirFes.O. 
(0.24% =0.8) has been carried out at 9300 Mc/s to study their ferromagnetic 
crystalline anisotropy and g-value. For all compositions a monotonic increase 
of Ki/M was observed as the temperature decreased, similarly to the curves 
obtained from the torque method. As a whole, however, considerably smaller 
values were obtained than the latter. g-value, also, increased as the tempera- 
ture was decreased for all compositions and showed larger variations as « was 
smaller. Qualitative explanation was tried from the view point of Clogston’s 
thermodynamic theory. 


§1. Introduction 


Recently, we have investigated Fe-Mn mixed ferrites (Mn.Fes-2O,, 
0.2<a%<1.15) and found that their ferromagnetic crystalline anisotropy has 
the maximum value in the region 0.6<a#<0.8'. A theoretical interpretation 
for this phenomenon on the basis of the one-ion model has been reported”. 

In order to study this behaviour deeply, the ternary Mn-Fe-Ni ferrite system 
has been investigated by Miyata in our laboratory. His experimental results 
and discussions on the ferromagnetic crystalline anisotropy which was obtained 
by the torque method have been reported®), and upon the same specimens as 
mentioned above we have experimented on the ferromagnetic resonance and 


tried to compare our results with the static data. 


§2. Preparation of Specimens 


(a) Crystals 

Single crystals were grown in air by Miyata using Bridgman method. 
Their compositions are «=0.20, 0.35, 0.52 and 0.75 in the expression with 
Ni,Fes_,O:. These values are the results from the chemical analysis. 
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(b) Polishing 
In order to polish the specimens spherically a device developed by Bond?’ 
was used. The grinder is shown in Fig. 1. 


emery paper 
“\ 


nozzle (05 $) flange 


Fig. 1. Grinder. 


An emery paper was affixed to the inner wall of the circular hole, in 
which specimens were blown with the compressed air through a nozzle. After 
grinding, the specimen was etched for a few minutes in the hot HCl solution 
of about 9 normal to remove the distorted layer of the surface. 


(c) Determination of the crystal orientation 

The orientation of a specimen can be determined accurately by the X-ray 
diffraction method. In order to rotate the crystal about a [110]-axis, the spots 
of the reflection from the (400), (440), (444) planes in the Debye-Scherrer 
photograph were pursued carefully. Its accuracy is within +0.5°. 

Then the specimen was mounted on the rod in the cavity without changing 
the orientation, and fixed with the special dental cement. It was coated with 
a thin film of polystyrene or silicon resin. 


§3. Experimental Apparatus 


(a) Microwave components 

The block diagram of the whole apparatus is shown in Fig. 2. The cavity 
of transmission type (TEios mode) was used. (Fig. 3.) | 

When the absorption line is strong the detection is possible by a simple 
video type; that is, by observing the decrease of Q-curve with the oscilloscope. 
But the absorption intensity of a specimen which contains ferrous ions 
becomes so weak that a high sensitive detector is necessary. 

The magnetic field was modulated by the town 50c/s and the phase- 
sensitive detector similar to Schuster type®) was constructed. (Fig. 4.) 

As the field modulation was applied, the microwave frequency had to be 
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stabilized. The A.F.C. network shown in Fig. 5. was used here. The micro- 
wave was frequency-modulated with 455 Kc/s crystal oscillator and the trans- 
mitted power through the cavity was amplified with the tuned amplifier. 


== A.F.C 


Field modulator 


Fig. 2. Ferromagnetic resonance spectrometer. 


spece men 


thermocouple 


ig. 3. Cavity of trans- 

eee je Fig. 4. Phase-sensitive detector. 

The amplified signal was discriminated and fed back to the repellor of the 

klystron (2K25) through the d.c. amplifier. The fluctuation from the centre- 
frequency was about +100 Kc/s. . 

The heterodyne frequency-meter is shown in Fig. 6. This is convenient 
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to observe the change of the resonance frequency of the cavity when the 
temperature is varied. Fig. 10 shows such a measurement. 


455 Ks 


oscillator 


Wyseon = Cavity 


Crystal 455KY/s 
detector amp. 


Fig. 5. A.F.C. network. 


Mixer 
(GAKS) 
fi 


U.H. F. =e k; 
oscillator 


R= -f, tf (n=23-26) 
zero beat F, = 4m (m+85~101) 


[| 


[FoF 
Crystal oe 


iver 
mixer Receive 


1Fi-4m] 


4Me/s 
drystal 
oscillator 


TE Cavity 
Directional 
Couple rm 


Fig 6. Heterodyne frequency-meter. 


9" 


656 


frequency Aange 


oe ae ae 340~ 404 "HC 


ST 


10K GLR) 


Fig. 7. U.H-F. oscillator. 


ae 
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plunger 


microwave 


U.H.F all. Coupling Loop (UNE input ) 
Fig. 8. Crystal mixer. Fig. 9. Mixer (6AK5). 
Meg 
9285 A 
> 
& 
290 8 
a 
5 
o 9215 © 
oO 
9270 
9265, 
too 150 ‘200 260 300 
T’K 


Fig. 10. Resonance frequency of the cavity vs. temperature. 


(b) Measurement of the magnetic field 

It is important to measure the magnetic field accurately. We used the 
fluxmeter with a search coil which was calibrated carefully, and also the field- 
meter of the proton resonance. ‘The electron spin resonance of D.P.P.H.® was 
observed frequently as a marker. 

For the field-meter of the proton resonance, the same circuit as Knoebel 
and Hahn’s”) was made. The flexible probe was convenient for our experiment 
and their frequency modulation method was not necessary in our case, where 
the magnetic field modulation was applied. The modulation amplitude of mag- 
netic field was varied from -++5 to +30 oersteds according to the absorption 


line width. 


§4. Analysis 


The resonance condition which has been derived by Kittel® is 
o=7 ett ? ( 1 ) 


where wo is the frequency in radian/sec, 7 is the gyro-magnetic ratio (ge/2mc), 
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g is the spectroscopic splitting factor, and 
Het=Hit+fil0)K1/M+fo(0)K2/M , (2) 
where @ is the angle of rotation in a crystal plane.” If we let 0 be the angle 


in the (011) plane between H; and the [100] direction, the following formulas 
are obtained, 


fi(0) =2—(5/2)sin? 0 —(15/8)sin? (20) , \ (3) 
f2(0) =sin? 0(1—6 sin? 0+(21/4)sin‘ 6) . 


Table. I. fi(@) and f2(@) in [100], [110] and [111] crystal axis. 


[100] [111] | [110] 

8 0° 54°44 90° 
fi (0) 2 —4/3 i172 
fo (6) 0 —4/9 1/4 


From Equations (2) and (3), He, Ki/M, and K2/M were determined with 
the measured values of the external applied field H. as a function of 6. Their 
temperature dependence was obtained from the data measured in the [100], 
[110] and [111] directions.* Also g-values at each temperature is given from 
Equation (1), or from the simplified formula: 

us f(Mc/s) 
Hes(oersteds) < 1.400 ’ 


where f is microwave frequency itself. 


(4) 


§5. Results 


The variations of the applied field required for resonance with crystal 
direction in the (011) plane are shown in Figs. 11, 12, 13 and 14 for four speci- 
mens respectively. The temperature dependence of K:/M, K2/M, and g-value 
for each compositions is shown in Figs. 15, 16 and 17. 

Two experiments on the single crystal of the Fe-Ni ferrite system have 
been reported by other investigators: D. W. Healy", and W. A. Yager et al!”. 
One of the two crystals which were studied by Yager et al is nearly equal in 
composition to our Nio.sFes 20:. But as to the temperature dependence of 


Ki/M, their results differ substantially from ours, on the other hand Healy’s 

* f2(@) takes zero at three values of @; one is, of course, @=0 and the other are real 
solutions of the next equation: 

1—6 sin? 6+(21/4)sin‘ @=0 . 

Solutions are given with @=26°44’ and 75°50’, so the data in these angles seem to simplify 
the calculation. But according to Artman’s consideration), in the neighbourhood of such 
angles the non-line up effect of the magnetization varies the resonance field considerably. 
So measurements in these directions are not necessarily available. 
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Ni 0:75 He 2.25 0, (O49 % 9) 


(Oe) 

200 A: 71'K 
CB , "kK 
$ 0: 290 
て 100 
: 
<= 

0 

-loo 

~200 

-300 

- 400 


Fig. 11. (Hres—Herr) vs. angle between M and [100] direc- 
tion in (011) plane for Nip 5Fes.25O4. 


Ni ost Fe 248 0; (047 %9) 


+300 Seats 
(Ge) Bi 19°K 
(as O: 290°K 
> 42006 
5 
< 
に 
= 
~ +100 
0 
-100 
-200 


Fig. 12. (Hres—Her) vs. angle between M and [100] direc- 
tion in (011) plane for Nio 52F ee. 4804. 


(0). NioasFe 2.65 0, (0.44 hp) 
A: 77K 
we O: 290°k 
= 
zx 
1 +100 
£ 
WwW 
| -100. 
-200 
301 
Fig. 13. (Hres—Hen) vs. angle between M and [100] direc- 


tion in (011) plane for Nio.ssFesz.65O4. 
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dog : Niigoore 2800, (0.50%?) 


(oe) 
4: 77°K 
©: 290°K 


Fig. 14. (Hres—Hen) vs. angle between M and ]100] direc- 
tion in (011) plane for Nip 20Fe2.s0Ox. 


Nk ageFe,. Bo 0, 


Nigar Feocs 0, | 


Nios2Fes.8 0, —>* 


se 


-s 


Moras Fe 225 0, 


~Ioo 


Fig. 15. Ki/M vs. temperature. 


NioroFe 2900i% 


Nio.as FeggsOu-* 2 


-200 


Nios2 Pe 2.4904 


Né 0.95 Fe225 Oy ——> 
=600 


Fig. 16. Ks/M vs. temperature. 
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23 


Nlo.75 Fe2.26 On 
SS 


22 Nio.as Fe265 Og 


ae 
Nio.s2Fes.48 0g 


Néo20 Fe 2.90 On 


21 


2.0 


Fig. 17. g-value vs. temperature. 


data are rather similar to ours, in spite of the different composition. The 
discrepancy may be attributed to the different procedure for growing single 
crystals. In our case, it is confirmed that the increase of Ki/M with the 
temperature decrease as shown in Fig. 15 is not caused by the cobalt ions?) 
included as impurities. 

The accurate measurement of the line width was not carried out, but at 
room temperature it was about 200 oersteds for Nio 1sFez 2504 and 300 oersteds 
for Nio,20Fe2 ss04 respectively. These widths became less than 100 oersteds at 
liquid nitrogen temperature. 

Figs. 18, 19 and 20 show how the Ki/M, K:/M and g-value vary with crystal 
composition. 


+200 


oO: 977°K 
A: /90°K 
O:290°K 


200 


Fig. 18. K.i/M vs. composition. 


18 


S. UsAMI and Z. FUNATOGAWA 


SaaS 9 T3270'K) 


Fe30, x NiFesOy 


Fig. 19. Ks2/M vs. composition. 


A 77°K 
0 290° 
2.4 


23 


2.2 


24 
20 
0 OF day 24 04 05 06 O71 og 0F 1.0 
Fes0, N Fes Og 
ーーー> エエ 
Fig. 20. g-value vs. composition. 
0.4 
O:authors 
«3 
re @ Healy + Johnson 
02 
or 


O Of 02 03 04 05 06 OF O88 09 ro 


I Nche,0, 


Fig. 21. g(77°K)—g(290°K) vs. composition. 
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§6. Discussions 


The static measurement by means of the torque-meter gives the similar 
tendency as our resonance data. Both results at room temperature and at liquid 
oxygen temperature are tabulated in Table I]. The saturation magnetization 
M was measured by Miyata and Nakamura in our laboratory for the same 
specimens. Fig. 22 shows the curve of 4Kivs. composition, where 4K, is the 


Table II. Magnetic anisotropy constant Ki, Kz in erg/cc, and anisotropy field 
Ki/M, K2/M in oersteds, as dependent on frequency and temperature. 


2. |temp.| K/M| Mm |K,(res.)| Ki(static)| 4K, | Ke/M |K,(res.)|K,(static) 
a9 | 2% | 413e°{ 460 | 46.3 | +12.3 | —6.0 | —100 hs 9 
(290°K) |(—147) | (440) |(—6.5) | (—6.7) |(+0.2) |(—100) | (—4) 
Ok le tee aio || 45.1N| 47.8. 4.8 | teol) —10 | —10, 
a (290°K) |(—134) | (878) |(—5.1) | (—3.3) |(—1.8) |(4100) | (44) _ = 
OR | 415, |..870> | 40.56] 44.4 | 3.8 | 270 -|.+10.| 12" 
a (290°K) |(—94) | (338) |(—3.2) | (—3.2) | © |(—103) | 3 | 5 
DT Spore ola-ds weheagenb=705ed) 223) bo=19 
© | aoe) |(—120) | (298) |(—3.6) | (—3.8) |(-+0.2) |(—120) | (—4) 


O: Ke (resonance) 
0 :Kal state ) 


A: BK, 


AK, (et 290°) 
x 


=) 


Fig. 22. Ki(res), Ki(static) and 4K, vs. composition at 
room temperature and at liquid oxygen temperature, 
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difference of the crystalline anisotropy between these two data: 
4K.:=K(resonance)— Ki (static) . 


Bozorth et al’ pointed out that in a Ni-Fe ferrite (NioFe216O.) the 
frequency dependence of the ferromagnetic crystalline anisotropy was observed. 
In our case, also, the same character has been observed in all compositions 


as seen in Table II or in Fig. 22. 


300, 


Nio-s2 Feoas 0, 


Nioao Fe2s004 
“8 
Fig. VSS Ak, vs. temperature for Nio 20F 2.3004 and Nio ssFes 4304. 


Bozorth et al regarded that this phenomenon was due to the relaxation 
associated with the rearrangement of electrons between ferrous and ferric ions 
minimizing the free energy of the crystal. 

From the above point of view, Clogston,'® also, tried to explain the tem- 
perature dependence of the line width, the ferromagnetic anisotropy and g-value 
in Nio.1sF' ez.2504, which were reported by Yager et al”), on the basis of thermo- 
dynamics. 

The temperature dependence of g-value in our study, too, may be accounted 
for by the thermodynamic theory. As for tne anisotropy, 4K: must be adopted 
rather than Ki(resonance) itself, since Ki(static) is not independent of the 
temperature. The temperature dependence of 4K: shows the same direction 
which the theory has expected. It is noticeable that at liquid oxygen tem- 
perature 4K: approximately depends linearly upon the composition (a#) and 
may be extrapolated to zero at the point of «=1.0. 

The Fe-Mn ferrite system which was investigated in our laboratory”.®).1% 
did not show such a frequency dependence as seen in the Fe-Ni ferrite system. 
At present the satisfactory explanation is not given for this difference. 


Acknowledgement 


The authors wish to express their gratitude to Dr. Miyata for specimens, 
static data and many useful discussions. They were helped also by many 
other people in our laboratory. This work was partly supported by the Scientific 
Research Expenditure of the Ministry of Education. 


15) 


Ferromagnetic Resonance in Ni,Fe3_,O, PAL 


References 


Z. FUNATOGAWA, N. MIYATA, and S. USAMI: J. Phys. Soc. Japan 14 (1959) 1583. 
J. C. SLONCZEWSKI: J. Appl. Phys. 32 (1961) 253S. 

N. MIYATA: J. Phys. Soc. Japan 16 (1961) 1291. 

W. L. BOND: Rev. Sci. Instr. 22 (1951) 344. 

N. A. SCHUSTER: Rev. Sci. Instr. 22 (1951) 254. 


A. N. HOLDEN, W. A. YAGER, and F. R. MERRITT: J. Chem. Phys. 19 (1951) 1319, 
H. W. KNOEBEL, and E. L. HAHN: Rev. Sci. Instr. 22 (1951) 904. 

C. KITTEL: Phys. Rev. 73 (1948) 155. 

J. R. BICKFORD: Phys. Rev. 78 (1950) 449. 

J. O. ARTMAN: Proc. Inst. Radio Eng. 44 (1956) 1284; Phys. Rev. 105 (1957) 62. 
D. W. HEALY: Phys. Rev. 86 (1952) 1009. 

W. A. YAGER, J. K. GALT, and F. R. MERRITT: Phys. Rev. 99 (1955) 1203. 

D. W. HEALY, and R. A. JOHNSON: Phys. Rev. 104 (1957) 634. 

R. M. BozortuH, B. B. CETLIN, J. K. GALT, F. R. MERRITT, and W. A. YAGER: 


Phys. Rev. 99 (1955) 1898. 
R. F. PENOYER, and M. W. SHAFER: J. Appl. Phys. 30 (1959) 315S. 


16) A.M. CLoGSTON: Bell Syst. Tech. J. 34 (1955) 739. 


H Sy e 7 
Co a 
wt 
1 ¥ ‘ 
s yf 8 
{3 x asf <4 
: ve OT ae 
bs au : RF A EE 
1 ices a 4 A W Vf 
Att mat Uf fi (IAG W 1 Ey : 
> wy Wt 4) A TAOIOG i eae 3 
Ei NY & “i weit 
Me, LL, IANS .W A bog eon 28 


ES vee! eT yr tee serene Wee a 
‘4 ’ 


os PTT We ts Ie ctr Gd 


wf sled hc ree turing aay <a 


| ve pont 2! few Kg le (eet an expinin oa pe 
ie oe Te re eh, thee fey dey ities te ding su 

の hap Pa by Upper 6 oy) ey tis tone of 
Syed ct = 5 
Thee Ursus tiks a os ¥ ie 

Tet bap ate Thee veal gaan, t Aa tat te ie 

iether ey Het scone) fn Hen Parcs Bae, 


rhe, PR is 
‘ SM notin 2h atl ha 


a 


e, ae 


[Science Reports of the Yokohama National University, Sec. I, No. 8, 1961] 


A Torsion Pendulum Magnetometer used 
for Thin Films 


By Yasuo GONDO and Hiroshi Konno 
(Received June 30, 1961) 


Synopsis 


In order to measure a small magnetic moment of thin films, a torsion 
pendulum magnetometer having a high sensibility is first constructed. Details 
of the principle and the construction of the apparatus are given. The magnetiz- 
ing field, in the range from zero to 10,000 Oe, applies parallel to the film plane, 
so that all specimens are magnetized to the saturation. The thickness and 
the temperature dependence of the magnetization in thin films have been 
determined by this method. 


§1. Introduction 


In the case of thin films, magnetization measurements are fatefully difficult, 
since, owing to their small volume, the magnetic flux or force is too small to 
detect. Several investigators have measured the saturation magnetization of 
thin films by various methods. These methods used for thin films are given 
in Table I. In both methods using the ferromagnetic resonance® and the 
Hall effect,» for example, it is necessary to saturate the film specimen 
perpendicularly to the film plane. Though it is usually assumed the demagnetiz- 
ing field is 4rM, but the effective one is to be very different due to the 
magnetoelastic anisotropy in the specimens with the large magnetostriction 
constant, which is caused from the remarkable streess anisotropy existing in 
thin films. If the stress anisotropy of the film is not estimated, therefore, 
the determination of the magnetization by those methods is rather unreliable. 

There are two inconsistent experimental observations on the thickness 
dependence of the magnetization in thin films. One says that it decreases 
with films thickness®) and the other insists that it is the same as a bulk speci- 
men®). Where these differences between them come from, is a question. One 
answer is that such a disagreement may be caused from the difference in the 
condition of the film preparation. While Crittenden and Hoffman* used a 
vacuum of the order of 10-*mm Hg, Neugebauer® used an ultra high vacuum 
of the order of 10-°mmHg or less. Another answer, however, is that the 
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Table I. Methods of magnetization measurements used for thin films. 


Methods Authors Remarks 
Direct method 
Induction \ * 
i) Ballistic galvanometer Sorensen,!) Reimer?) | <800 Oe. Several Films must 
be used together. 
ii) Loop tracer Crittenden-Hoffman®) | <400 Oe. Signals are amplified. 
Body force : 
i) Magnetic balance Jensen-Nielsen?) <10,000 Oe. Non-uniform field. 


Demagnetizing field ~4nM. 
<10,000 Oe. Non-uniform field. 


ii) Torsion pendulum mag- | Gond6-Konno 


netometer (this report) Demagnetizing field is zero. 
iii) Torque magnetometer | Neugebauer?) <10,000 Oe. Uniform field. 
Angle between film plane and 
field ~7/4. 
Ferromagnetic resonance Seavey-Tannenwald® | Resonance field fixed. Demag- 


netizing field is ~4nM. 


Indirect method 
Galvanomagnetic effect 


i) AR effect Colombani,”) ~10,000 Oe. Curie point. 
Kuwahara®) 
ii) Hall effect Coren-Juretchke®) ~10,000 Oe. Demagnetizing 
field ~4nM. 
Magneto-optic effect 
i 10) 7 2 . 
i) i ee To obtain magnetization curve. 


difference between their methods of the magnetization measurement may be 
an origin of the disagreement. The former, the induction method, used a 
magnetic field of less than 400 Oe, which was too weak to saturate the speci- 
mens, and the latter, the torque method, used up to 10,000 Oe which was 
enough to saturate them. 

In order to make clear whether these differences are due to the measuring 
methods used by them or not, it has been desirable that the films are prepared 
in a vacuum of 10-6mm Hg and measured by a more improved method. So 
we have constructed a modified torsion pendulum magnetometer, one of the 
body force method"), to determine the magnetic moment of thin films and 
measured field- and temperature-dependence of the magnetic moment, too. 


§2. Principle of Measurement 


A thin film is placed at the middle point O of an air gap between circularly 
capped pole pieces of an electromagnet. Here the normal of the film plane 
is perpendicular to the axis of a set of the pole pieces, the x axis, and the 
specimen oscillates transversely to the axis, keeping the direction of the normal 
along the y axis. This is shown schematically in Fig. 1. Because of the shape 
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anisotropy, the magnetic moment of the film should be in the film plane along 
the x axis. Therefore, when the field applies parallel to the film plane, the 
film is easily magnetized to the satura- 
tion. It can be considered the magnetic 
poles on the spherical pole piece are 
converged together in the center of the 
sphere. It is obvious that the field at 
the point (0, y) on the y axis is given by 


2D 
A.=m ーー 
(D?-+-y?)3/2 (1) Fig. 1. The illustration of the mag- 
and netic restoring force F exerted on 
H,=0, the film in the field. 


where m is the magnetic pole strength and D is the distance between O and 
the point for the poles to be converged together as shown in Fig. 1. Since 
the x component of the field at O, Ho, is 2m/D?, the above equation is written as 


If y/D<1, eq. (2) gives 


=H {i— + (#)} (3) 
= ca 
When the film is displaced from the middle point O along the y axis, the 
magnetic force exerted on the film at the point (0, y) is given by 


0H. __ 3 


ay pe ey » (4) 


F,= Mv 


where M is the magnetization and v is the volume of the film. This magnetic 
force is always opposite to and proportional to the displacement of the film 
from O, so that, when its motion is restricted only along the y axis, the film 
makes a simple harmonic oscillation along the axis. 

When no magnetic field is applied, the period of a pendulum system is 


given by 


seed oe ke ei (5) 


where J is the moment of inertia and ko is the restoring torque constant of 
the system. When the magnetic field is then applied, the ordinary elastic 
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restoring force is superposed with the magnetic restoring force, as given by 
(4), and consequently the period of the system changes remarkably. Then 


we have the period 7 as 


= I 
T= y/ >t, (6) 


where km is the torque constant due to the magnetic force. Here it is given 
by 


kn=F,?2 = —— = (7) 
where / is the radius of gyration of the system. From (5) and (6) we have 
he = Rm 
fete Mes FA 8 
les Ro (8) 


This relation is shown in Fig. 2. Since the magnetic moment Mv is very 
small, km is not large for thin films. 
Hence, in the case of measurements for 


thin films, it is necessary that we make 
eS ky) be small as much as possible, to 
Os obtain a high sensibility. Consequently 
the torsion type should be adopted for 

? the pendulum system. 
o 4 2 3 2% Finally the following equation is 


Fig. 2. The relation between the mag- derived from (8) 
netic restoring constant km and the 


period T. AMoHy=(=) 1 : (9) 


where 


ee 3 Bt ST at 
Dko 422 1D? 


A (10) 


is a constant determined by the geometrical condition of the apparatus. So 
the intensity of the magnetization of the specimen can be estimated from eq. 
(9), as a function of the field strength Ho, when A is known for the apparatus 
and 7, T and Hy are measured, respectively. In this report afterwords HA 
will be written as H in abbreviation. 


§3. Preparation of Specimens and Measuring Apparatus 


The circular films of nickel and other metals are prepared by evaporation 
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from a tungsten filament on to cover glasses for microscope held at temperature 
of 100 to 200°C. The diameter of the films is 0.8cm. The vacuum condition 
during the evaporation is the order of 10->mm Hg. They are coated with 
SiO immediately after the evaporation to prevent oxidation. 

The apparatus used for the 
measurement is essentially a torsion 
pendulum and the general assembly 
is shown in Fig. 3. 

The pole piece of the magnet 
is 6cm in diameter and 3.6cm in 
radius of the spherical plane. The 
field distribution along the y axis 
always decreases parabolically with 


y near the center as shown in Fig. 
4, so in this region the condition Fig. 3. Schematic diagram of the apparatus. 
V: To the vacuum pump. W: Torsion 
wire. B: Counterbalancer. M: Mirror. 
complete. A: Arm. D: Dewar flask. F: Film speci- 
men. L: Lamp. P: Photo-tube. N, S: 
Pole pieces of the mganet. 


for simple harmonic oscillation is 


The torsion wire of the pendu- 
lum is made of a constantan wire 
of 604 in diameter and 10cm in length. It is desired that the arm of the 
pendulum is long enough, because the specimen has to do a linear motion 
approximately in the range of the amplitude. In this case its length is 20cm 
and the amplitude is limitted less than 0.5cm. This pendulum arm consists 
ay of a thin aluminum pipe hung 

horizontally by the torsion wire and 
a paper pipe and a specimen holder 
connected vertically with the above. 
The specimen is fixed on the holder 


by high vacuum silicon grease and 
Fig. 4. The relative field distribution along settled at the center of the air gap. 

Soc et tat There is a counterbalancer on the 
other side of the arm. The arm length from the torque wire to the specimen 
is 20cm and then km is 2.5x10-2 dyne-cm. This pendulum system must be 
placed in a vacuum in order to avoid the disturbance of air. To examine the 
effect of temperature on magnetization, a Dewar flask is prepared as shown 


in Fig. 3. 
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§4. Measuring method and Results 


When the pendulum system oscillates, a transit of the light beam reflected 
by a mirror, attached on the head of the pendulum arm, is transformed into 
an electric pulse signal by a photo-tube connected with a recorder. Comparing 
the record of a sequence of the transits with the standard time signal from 
JJY*, we can measure the period of the system. 

We set the specimen on the 
pendulum and make it oscillate, 
at first applying a strong field. The 
period T with the field is measured, 
then its change with decreasing the 
field is examined. Examples of the 
observed results are shown in Fig. 5. 
The period becomes longer with 
Fig. 5. The changes of the period with the decreasing the field, because of re- 

field. ducing the magnetic force exerted 


on the specimen. The period changes in the range from a few seconds to 50~60 
seconds. This is measured with an error of less than 0.5%. Provided the period 
To, without field, is also measured, using (9) we obtain the value of the magnetiza- 
tion M as a function of H. The typical ones are shown in Fig. 6. Thus the 
hysteresis curve of the specimen as 
well as the saturation magnetization 
is determined accurately. It is 
noted that we need the field of 
2,000Oe or more to saturate the 
film. If the applied field is not so 


strong, therefore, we can not deter- et SS ey Ve 
mine the saturation magnetization. 4 Koe 
Fig. 6. The magnetization curves obtained 


Let us consider the effect of F 
from Fig. 5. 

para- and diamagnetic bodies in the 
apparatus. The magnetization, observed by induction or body force method, is 
the total magnetization of the specimen itself, the substrate and the specimen 
holder. In the case of thin films, the magnetic force due to the ferromagnetic 
body is not so large that we can not neglect the force due to the para- or 
diamagnetic substrates. To obtain the exact results, therefore, we have to 
subtract them from the observed values. If the diamagnetic contribution is 


larger than the ferromagnetic one, the pendulum goes away from the center and 


* The transmitting station of standard frequency in Tokyo. 
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does not make oscillation. It should be also considered that the condition of 
oscillation changes remarkably owing to the effect of the electrostatic Coulomb 
force between the charges on the arm and the wall of the vacuum system. 

When the air gap of the magnet is 20 mm, the field H is provided up to 
about 10,000 Oe, which is usually enough to saturate the magnetization of thin 
films. The determination of the field H is made by measuring the intensity of the 
exciting current with a great care, because there is a hysteresis in the magnet. 
The determination of the field, however, can be done within the error of 1% 
or less, even in the region of lower fields, if the same descending hysteresis is 
used. 

This torsion pendulum magnetometer is enough to observe the saturation 
magnetization of thin films, their magnetization process and its temperature 
dependence, too. The apparatus has the sensibility, capable of the measure- 
ments down to 10A in thickness. The experimental error for the observation 
are 10% or less in the case of saturation magnetization measurements of 
nickel films of 20A in thickness. The error is almost due to the one of the 
estimation of the film thickness, that is the volume, and the error of the 
obtained value of the magnetic moment itself is less than a few per cent. 
By this method it has been found that nickel films have large anisotropy 
perpendicular to the film plane which may be caused from the tensile stress 
in the film plane, so they are scarcely saturated with the field of less than 
3,000 Oe; and that the saturation magnetization is the same as bulk nickel 
down to a thickness of 25A; and that in the very thin films of about 20A 
or less in thickness the saturation magnetization and its temperature depend- 
ence seem to be due to the superparamagnetism. These investigations are now 


in progress.!®) 
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Lichtenberg Figures Observed on Color Negative 
Films in Low Atmospheric Pressure 


By 
Bunji ARAI 


(Received June 10, 1961) 


Synopsis 


The experimental studies were performed as a part of the serial study to 
elucidate structure of Lichtenberg figures and the process of their formation. 
In the experiments we used the color negative films were used instead of 
the reversal color hitherto used. 

The Lichtenberg figures on the usual black and white dry plates under 
various pressures were investigated by J. Stekolnikof and K. Riaschenzef and 
others. 

The color negative was newly used because of its sensitivity (ASA 32° 
or 50°) superior to the reversal color (ASA 12°) hitherto used with the view 
of finding important phenomena otherwise overlooked. We not only observed 
color Lichtenberg figures produced by the impulse of definite wave form and 
crest voltage under various atmospheric pressures. 

Observation was made mainly under the reduced pressures and the variation 
due to pressure change, of the size, form and color of the figures was investi- 
gated. 

The size of the positive and negative figures grew as the pressure was 
reduced, and the color of streamers showed a tendency to grow thick and bluish 
violet tint (yellowish orange on color negative films). While in the gathering 
fan-shaped negative figures we found the same tendency, the rate of their 
growth due to pressure decrease was more remarkable than the positive, and 
their red streamers on the periphery (cyan-blue on the color negative) became 
more remarkable as the pressure was reduced. The red streamers which 
appeared only very slightly under a normal pressure became very distinct 
under reduced pressure. Besides, inferring from the multiple layer construction 
of color film and the color observed on the figures, it might be concluded that 
positive figures developed on the top layer of the emulsion only, while 
negative figures had such a structure that the prong of their streamers crept 
as far as to the bottom emulsion layer. 

As for the rate of their growth, the negative figures showed superiority 
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to the positive when the pressures was remarkably reduced. Namely, we found 
a tendency of the ratio to approach from R+/R->2 under normal pressure to 
1 or Rt R- as the pressure was reduced. We confirmed the fact that the 
effects of pressure reduction gave dominant effects on the extension of 
streamers of both figures. 

The tendency of J. Stekolnikof and K. Riaschenzef’s curves was different 
from that of the writer’s. 

The study on the color Lichtenberg figures under further decreased pressure 
was made to trace the relation among the pressures, the creat voltages and 
the variation of both figures. In these cases, the color of the streamers of 
both figures was very weak bluish violet, and we could not find any trace of 
red streamers on the tips of the negative figures. The critical pressure was 
estimated to be about 7cm Hg, under 5.5kV of the crest voltage supplied. 
From these facts we might conclude that both Lichtenberg figures were formed 
in such a case only on the surface of emulsion layer. 


§1. Introduction 


The principle of ‘‘ New instrument....New result’’ is the dominant theme 
of this work. The writer tried to obtain Lichtenberg figures by applying 
impulse voltages to color films instead of the photo-sensitive materials previ- 
ously used. These color films were used with the new instrument. The 
writer previously reported about positive and negative Lichtenberg figures 
which were obtained by supplying an impulse voltage to these color films 
(multiple layered, reversal type) under atmospheric pressure’. 

The reason for using color negative films in those experiments, was that 
it had greater sensitivity (ASA 32° or 50°) than the usual color reversal films 
hitherto used (ASA 12°) and, it had multiple layers. From this point of view 
it was hoped that more important phenomena could be investigated that had 
been investigated previously. 

On the usual black and white dry plates, the studies of Lichtenberg figures 
under the changing atmospheric pressure were made by J. Stekolnikof and K. 
Riaschenzef”) and others*®). In this work we investigated Lichtenberg figures 
produced on color film under reduced pressures. 


§2. Experimental apparatus and Procedures 


We took up the method which C. E. Magnusson“) had used for obtaining 
both positive and negative Lichtenberg figures by the same crest impulse 
voltage at the same time. 
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(1) The Klydonograph camera 


The special camera was one which the author had made. This camera was 
made of a bakelite plate and a round bar, a brass plate and bars, glass plates, 


pieces of pasteboard, vis and nuts etc. The details of the camera are shown 
in Fig. 1. 


Fig. 1. Details of special camera designed for this experiment. 

C : Brass plate (8.3 cmx6.0 cm x0.15 cm). 

G, : Clean glass plates (8.3 cmx6.0 cmx0.13 cm, gelation removed dry 
plates. These glass plates were removed in the preliminary ex- 
periment when dry plates were used instead of films). 

F : Fuji Color negative cut films (a quarter of cabinet size), Fuji Al dry 
plates in the preliminary experiment. 

P, Q: Brass bar electrodes (0.30cm? in diameter x3.0 cm long), screwed 
in the bakelite plates. 

Bs : Bakelite plate frame (0.30 cm thick). 

Rp, Rg : Bakelite round bar lock nuts (1.4 cm®%). 

M, N: Terminals for high tension and earth side. 

X : Double cover of mercury lamp (pasteboard ok 

Y, Z: Pasteboard side frame for camera. 

Rp: Rubber band. 


The metal plate (brass plates, 0.15cm thick) was put in the center of the 
camera, and glass plates (clean glass plates after removing the emulsion from 
dry plates, 0.13cm thick) were put at either side of this metal plate. The 
color negative films were put outside the glass plates. The films were cut out 
in quarter size from one sheet of cabinet size color negative cut film. Two 
sheets of quarter size film were tested at a time in the experiment. The 
celluloid bases of the films were placed in contact with the metal plate, and 
the emulsion surfaces in contact with the bar electrodes outside this camera. 
The bar electrodes were screwed into the bakelite plates forming both side 
of the frames of this camera. Reference to Fig. 1 again. This camera was 
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made to have perfect electrical insulation, and to be light tight. The metal 
plate put in the center was intended to have unequally electric fields on its 
both sides between the plate and the bar electrodes with the color negative 
films inserted. 

The contact pressure of the bar electrodes was usually kept constant, by 
pressing the outer cover and the bakelite plate frames together with two 
rubber bands. 

(2) Equipment for changing atmospheric pressure 

The special camera loaded with the color negative films was put in a 
vacuum desiccator with three holes in it. The air pressure in this desiccator 
was lowered extraction through the trap by the rotary oil pump. The pressure 
was read by the Bourdon gauge, and the open and closed tube mercury 
manometers. The pressure was adjusted by the leak cock. Refer to Fig. 2. 
(3) High tension impulse source and circuit 

The high tension electrical equipment was the same as the impulse circuit 
hitherto used. It was supplied from half of the Gaiffe-Gallot et Pilon’s D.C. 
high tension apparatus!).»), Here, the high tension condenser of the impulse 
circuit was 5yF (the dielectrics were injected with SHIBANORU, and the max. 
D.C. working voltage was 20kV). R andr were nearly equal to 1.11052, 
and 0.611032 respectively. Refer to Fig. 2 again. 

(4) Experimental procedures 

The spherical spark-gap Gi: (2.5 cm?, brass ball) was kept at the length of 
0.30cm with a few exceptions when it was adjusted to 0.20cm or 0.36cm. 
The calculated crest voltage Vn between P and Q corresponding to the 0.30cm 
spark-gap is 11.0kV (7.4kV or 13.0kV for the other gaps) at 25°C and 76.0 cm Hg. 
The other spherical gap Gz (2.0 cm?, brass ball) was kept at the length of 4.0 cm 
through the experiment. Both gaps were illuminated with the ultra violet 
rays from a quartz mercury vapour lamp to minimize fluctuations of the spark- 
ing discharges. The values Vm of crest voltages observed with the aid of a 
high tension voltmeter (for television use) were confirmed to be in good agree- 
ment with the calculated values Vn. 

The color negative films used for these experiments were cut films of 
cabinet size (UCN PC/135-02). These had been made in the Research Laboratory 
of the Fuji Photo-Film Co. Ltd., and were kindly presented to the author. 
The films were cut to a quarter size from the cabinet size film with a jig in 
a dark room, and loaded in the special camera. They were exposed with a 
constant crest voltage at the desired ambient pressure. After the impulse 
was supplied to the films, they were taken out of the vacuum desiccator, 
unloaded from the camera in the dark room, and put in the black paper of 
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the tin foil envelope. All the color developments and the color prints of these 
color negative films were submitted to the Research Laboratory of the Fuji 
Photo-Film Co. Ltd. for processing. 


Fig. 2. High tension impulse circuit to special camera and equipment for 
changing atmospheric pressures. 


T : High tension transformer. 

K,, Kz: Kenotrons (filament voltage 10V). 

C,, Cz: High tension condensers. 

G; : Sphere spark-gap (2.5 cm¢, brass ball). 

Gp : [2 CU ”). 

Cs : High tension condenser (5 pF). 

Rs : Lig. resistor (dilute solution of H2SO., 105.2 x 2). 


Rae ” ri ee [2 OR 
fend ” » for damper of impulse wave (solution of CuSQ,, 
0.61 kf). 


Vr: High tension voltmeter (for television). 


K Special klydonograph camera which was put in yacuum desiccator. 
a. S) serap: B: Bourdon gauge. 

O Open tube mercury manometer. 

Cc Closed 7 ” ” 

L Leak cock. 


Fuji Al dry plates were used for the preliminary experiments, and the 
tendency of variation with pressure was observed. In this case, the glass 
plates between the bar electrodes and the metal plate in the camera were 
removed. The dry plates of one lot were developed simultaneously by the 
author himself with the developer DK (20°C, 5 minutes), and the fixing bath 


FF-H, (20°C, 10 minutes). 
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§3. Results of Preliminary Experiments 


In the preliminary experiments, we observed the variation of Lichtenberg 
figures in various reduced pressures with the black and white dry plates (Fuji 
Al). The value of crest voltage supplied on the films or dry plates was limited 
by the area of the emulsion surface and the degree of reduced pressure. 

Photo. 1 show the reprint of a series of the typical results of the pre- 
liminary research tried with the Fuji Al dry plates at a constant crest voltage 
of 11.0kV. We understood that the size of positive dendritic figures show 
the enlargement, and each streamer was longer and wider in dependently as 
the atmospheric pressure decreases in spite of the same crest voltage supplied. 
In the negative figures, we observed the same tendency toward enlargement. 
But the rates of their enlargement for the reducing pressure were remarkable. 
In this, our results differ from those of J. Stekolnikof and K. Riaschenzef as 
will be discussed later. 

Fig. 3 shows the diagram of the relation between the decreasing atmos- 


pheric pressure and the radius of both figures (the length of each streamer) 
Fuji Al Dry plate (9105-459A). 


Radius of figures in cm RY = 
Fig. 3. Relation between the decreased Fig. 4. Relation between the reduced 
pressures and the max. radius of both pressures and the ratio of radii of 
figures in the preliminary experiments. both figures. The notations are 


the same as those in Fig. 3.. 


Fig. 4 shows the relation between the reducing pressure and the ratio of 
the radius of both figures (R+/R-). The parameter in this diagram was crest 
voltage. In each case, we found the very interesting fact that, as the pressure 
is reduced, the ratio of the length of each streamer R+/R- approaches 1 from 
R*/R- >2 at normal pressure. That is, R+ becomes nearly equal to R-. The 
extension of the streamers of the negative figures was greater than that of 
the positive figures as the pressure was reduced. 


Lichtenberg Figures Observed on Color Negative Films 37 


We confirmed the fact that the effect of reducing pressure greatly effected 
the extension of streamers of both positive and negative Lichtenberg figures. 


§ 4. Results and Discussion 


Photo. 2 shows the typical series of Lichtenberg figures reprinted on a 
black and white photographic paper (chloried or gaslight paper, because the 
color reprint was not ready at the time of this report), which were obtained 
by supplying a constant crest voltage (11.0 kV) at various reduced pressures on 
the color negative films. The enlargement and variation of the streamers on 
both figures has nearly the same tendency as the result of the preliminary 
experiments by the Fuji Al (orthochromatic) dry plates. 

The color of the streamers themselves was blue-violet in the positive 
figures irrespective of pressure, and the halo in contact with the neighbouring 
bar electrode was white-blue. In the negative figures, the color in the region 
of the neighbouring bar electrode showed white-blue, but the prong near the 
circumference in the figures was yellow-orange (the color of this part in the 
color negative film was cyan-blue). These yellow-orange or red streamers 
showed that the bottom emulsion layer (the nearest layer of the celluloid 
base) of this color film was excited by the negative streamers. 

Therefore, we understood that the negative streamers crept into the 
bottom emulsion layer from the surface. From characteristics of the multiple 
layered color film, these facts give an actual proof that the positive figure 
was formed at the emulsion surface, while the negative figures had a volume 
construction. This phenomena also appeared under normal atmospheric 
pressure”).*), but it appeared very remarkable with very low crest voltage 
under reduced pressure. 

The diagram of Fig. 5 shows the relation between the streamers of both 
figures in Photo. 2 at reduced pressures. They are nearly equal to the results 
of the preliminary experiments shown in Figs. 3 and 4. The diagram of Fig. 6 
shows the results of J. Stekolnikof and K. Riaschenzef’s observation. Their 
curves are different from that of the writer’s. 

Their curves are as follows: 

1) Plotted from four observed points, they were not curved in one 
direction, and had a point of inflexion in them for the positive figures. 

2) The rates of the enlargement of the negative figures were very 
small compared with our results. 

We think that their observations resulted from the following causes: 

1) Their equipment is different from our apparatus. Mainly, the arrange- 
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ments of the electrodes were different; their equipment could not be used to 

obtain both positive and negative Lichtenberg figures, at the same time, and 

the constants were different from those of the writer’s impulse circuit. See 
Fuji Color Nega. film Figs. 7 and 8. 

(UCN PC/135-02) 2) The photo-materials used by them in those 

Ratio both pees days were probably low and narrow in sensitivity 

compared with our present photo-materials. 


And also we think that their different results 


ee had probably been caused by the combined effects 
ae because of the preceding two reasons. 
5 “i As our results show in Figs. 3, 4 and 5, the 
a curves are smooth and have a single direction 
20 


with no point of inflexion. We believe our curves 
are correct, because we made more observations 


and obtained more data than Stekolnikoff and 
Radius of figures incm  Riaschenzef. 

Fig. 5. Diagram showing the Next we tried the following experiments in 
relation between the reduced 
pressures and the size of both 
figures, and that between the with the same crest voltages. In this case, the 
pressures and the ratio of the “glass plate and the color negative film inserted 
radii of both figures. x ; Min. 
radius where red streamers 
appeared. plate electrode were put on one side only. Refer 


which each figures was obtained independently 


between one side of the bar electrode and the 


MRE の 2 WD ME eS CPS! 6 RE TT 
Rinmm Rinimm 
Positive Figuren Negative Figuren 
Fig. 6. Diagrams showing the results of J. Stekolnikof & K. Riaschenzef’s 
observation. 


to Fig. 9. Then we observed Lichtenberg figures obtained after supplying the 
constant impulse voltage at various reduced pressures. In this case, the crest 
voltage was ajusted by the spherical gap G: to about half that of the former 
case, Photo, 3, 4 and the diagram of Fig. 10 show the results. In the 
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Fig. 7. J. Stekolnikof and K. Riaschenzef’s equipment. 


diagram, our curves show the same R-v'R 
tendency as that described above, which area Fe 
again is different from that in the report x 3 
of Stekolnikif and Riaschenzef.” 0 
The final decision as to the correct . 
results must be left to the next project. Fig. 8. Impulse circuit of the above. 


This will consist of studies on color Lichtenberg figures formed at still lower 
ambient pressure in order to trace the relations among the pressure, the crest 
voltage and the size of both figures. 


Photo. 5 and 6 were obtained when we decreased the ambient pressure to 
6.8cm of Hg. These are the photographs printed from the color negative 
films. In these cases the color of the streamers of both figures showed very 


Fig. 9. Arrangement to take 
photographs of positive and 
negative figures inde- 
pendently. In this case, the 
glass plate and the color 
negative film or the dry 
plate inserted between one 
of the bar electrodes and 
the brass plate electrode 
were removed. 


weak bluish violet, and we could not find the marks 
of red streamers on the tips on the negative 
figures. We observed how the critical pressure 
under which Lichtenberg figures changed in such 
a way varied with the crest voltage. We found 
that the critical pressure is about 7cm Hg, and 
11.0kV for the value of the crest voltage. 

Photo. 7 and 8 are the photographs of Lichten- 
berg figures obtained at 6.8cm of Hg and printed 
from the color negative films which we obtained 
by inserting the films in the inverse direction, 
that is, with the celluloid base touching the bar 
electrode, and the emulsion surface in contact 
with the side of the metal plate electrode. In this 
case the figures were all recorded as yellow-orange 
streamers. From these facts we deduced that 
both positive and negative Lichtenberg figures 
were formed only on the surface layer where the 
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bar electrodes contacted the films. 

The elongation of streamers and increase of brightness of the figure can 
be explained to be due to the increase of the mean free path of gas (air) 
molecules as well as to the growth of corona discharge under reduced pressure. 
We might infer that the difference 
in the form and elongation of 
streamers in positive and negative 
figures would show a difference in 
their formation mechanism and also 
in the nature of the charged particles 
contributing to the formation. 

We are now performing some 
preliminary experiments on the 


relations between the colored Lichten- 


Radius of figures in cm 


‘ ; het berg figures and the form of the 
Fig. 10. Diagram of the relation between 


the reduced pressures and the elongation 
streamers, the arrangement of Fig. 9 the wave form of the impulse are 
being used. 


impulse voltage. The figures and 


photographed by means of the hand 
made special camera and the high speed cathode ray oscilloscope (Synchroscope- 
SS-5301) at the same time. 

These results will give a dominant analytical solution to the construction 
of Lichtenberg figures and the process of their formation. 


§5. Conclusion 


Summarizing the experimental results: 

1) The pressure decrease of the atmosphere has a major influence on the 
size of positive and negative Lichtenberg figures; the elongation of streamers 
along the surface depend on the atmospheric pressure. 

2) The elongation of the negative streamer is more remarkable than of 
the positive figures under the same condition and we might infer that the 
change in the size of the streamer due to pressure decrease is continuous. 

3) The circumference of the negative figure creeps into the deepest layer 
of the emulsion and thus three dimension structure is formed while the positive 
figure is formed only on the surface of the emulsion. 

4) About the negative figure, we could confirm the circumferential 
development into the deeper layers even under the normal pressure when 
excited by a much higher crest impulse voltage. Under reduced pressure it 
appears in the deeper layers when a lower crest voltage is applied, 
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5) Below a certain critical pressure and with a given voltage, the three 
dimensional structure of the negative figure disappears, and both figures are 
abruptly formed on the emulsion surface only. 

From the above experimental facts, it is suggested that the low crest 
impulse voltages which could not be analyzed with the klydonograph previously 
used, can be detected as enlarged and bright figures by means of a special 
klydonograph under reduced pressure. 

A further project to make a still low lower pressure klydonograph is being 
undertaken. 

This report was prepared from lectures the author gave at the Special 
Meeting on Electric Discharge Physics, held by the Physical Society of Japan, 
in Hiroshima, on Oct. 9, 1959. 
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Fuji Al Dry Plates (9105-458 A) 
Positive figures Negative figures 
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Photo. 1. Series of both positive and negative Lichtenberg figures under various 
reduced pressures and a constant crest voltage (VO) 


Plate II B. Aral: 


Fuji Color Nega. films (UCN PC/135-02) 
Positive figures Negative figures 


p [cm Hg] 
76.1 
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28.1 
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Photo. 2. Series of both positive and negative Lichtenberg figures (contact print 
of Fuji Color Negative) under various reduced pressures and a constant crest 
voltage (Vn== Vm=11.0 kV). 
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Fuji Al Dry Plates (9105-589 A) 
Positive figures 


p [cm Hg] 


46,1 
28.1 


58.6 


48.1 


38.5 WM 


Photo. 3. Series of positive figures under various reduced pressures by the 
arrangement of Fig. 9. (Vn= V4 kV 
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Fuji Al Dry Plates (9105-589 A) 
Negative figure 


p [cm Hg] 
76.5 


68.7 


08.7 


48.3 


27.3 


20.0 


13.6 


Photo. 4. Series of negative figures under various reduced pressures by the 
arrangement of Fig. 9. (Vn Vm=6.0 kV). 


B. ARrai: Plate V 


Fuji Color Nega. film (UCN PC/135-02) 


Positive figures Negative figures 


Photo. 5. Printed from the Fuji Color Nega. films under a crest voltage 11.0kV 
and pressure further reduced (6.8 cm Hg). Color of both figures was blue. 


Photo. 6. Under the crest voltage of 13.0kV and with data the same as those 
in Photo. 5. Also color was blue. 


Plate VI B. ARAI: 


Fuji Color Nega. film (UCN PC/135-02) 
Positive figures Negative figures 


Photo. 7. Printed from the Fuji Color Nega. films under the crest voltage 11.0 
kV and the pressure further reduced (6.8 cm Hg), the Color Nega. films inserted 
in the inverse direction. Color of both figures was yellowish orange. 


Photo. 8. Under the crest voltage 13.0 kV with other data the same as those in 
Photo. 7. .Also color was yellowish orange. 
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Some Contribution to Lattice Model Theory 
of Polymer Solutions* 


By 
Tatsuro WATARI 


(Received April 20, 1961) 


A systematic method of constructing the combinatory factor 
for configurational partition function of an assembly on lattice is 
given in Chap. II. This method, called pseudo-assembly method, 
not only makes the construction of combinatory factors easier, but 
also clarifies their natures. 

In Chap. III, it is shown that the pseudo-assembly method is 
very efficacious when applied ‘to the construction of combinatory 
factors of chain polymer solutions, and that it gives configurational 
partition functions which lead to an improved lattice model theory 
of chain polymer solutions, giving the second osmotic virial coef- 
ficient in terms of shape factors introduced in order to specify the 
shape of a polymer chain. 


I. Introduction 


A. Cluster Variation Method 


In statistical thermodynamics, there are a large number of problems which 
can be reduced, after suitable rewritings, to those of assemblies on lattice. 
Although there remain, especially in the case of fluids, many problems in such 
rewritings themselves, once the problems have been reduced to those on 
lattice, one of the central problems is to find out the configurational partition 


function: 
Vi >Gs exp [—E;/kT ] ’ 


where G, represents the number of distinguishable configurations of an assembly 
of lattice elements on a given lattice when the interaction energy is Es, Rk is 
the Boltzmann constant, and T is the absolute temperature. Lattice elements 
represent, for example, atoms in the case of alloys, magnetic spins in the case 
of ferromagnets, molecules and holes in the case of pure liquids, and solute 


* Submitted to Tokyo University of Education to fulfil the requirements for the 
degree of Doctor of Science, under the title of Pseudo-assembly Method and its 
Application to Lattice Model Theory of Polymer Solutions, on Nov. 18, 1959. 
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and solvent molecules in the case of solutions. 
Usually, the total interaction energy is tacitly assumed to be the sum of 
those due to the nearest neighbors. In such a case, the partition function may 


be written as follows: 


Z= 2 GUN, {ysa, {hr}) exp [—E({Ni}, {ys.nt, @/RT] , (A.1) 
Jy 
where Nj; represents the number of lattice elements j, Yj.x the probability of 
finding a pair of cluster of lattice elements j and & on a pair of neighboring 
lattice points, G the number of distinguishable configurations of the assembly 
satisfying the set {y;,x}, and E is the total interaction energy of the assembly. 
Furthermore, in the above equation, h, represents the r-th index specifying 
the structure of the lattice. When the lattice with N lattice points has oN 
lattice bonds and the coordination number of the lattice is 2w, w may be assigned 
to the first specification index hi. It must be noted that, while E is determined 
by giving hi=o, G is determined uniquely only when the full set {h,} is given. 

In the case of a one- or two-dimensional lattice, the rigorous estimation 
of the partition function (A.1) has been carried out with the so-called matrix 
method by Onsager and others.'!) However, the three dimensional problem has 
been left unsolved by the similar method owing to the difficulties in carrying 
out the calculation practically. 

There has been another method more direct in its approach to this problem. 
Though this method has given only the approximate estimation of the partition 
function even in the case of a two-dimensional problem, it is easily applicable 
to the three-dimensional problem. The original type of this method owes much 
to the works of Bethe,®” Fowler,*) and Takagi.‘) In this method, i= is used 
as the sole specification index of lattice structure. Hence, G in Eq. (A.1) is 
constructed in terms of {N;}, {y;,.} and w. In the present paper, this method 
is called the first approximation method. 

Further developments have been made in this methcd by taking into 
account additional specification indices corresponding to the existence of higher 
correlations of lattice points. In the case of a triangular lattice with oN 
lattice bonds and vA triangular cells, for example, the existence of triangular 
cells in the lattice structure is taken into account by assigning v to the second 
specification index hs. Accordingly, the probability 2;.k.i of finding a triangular 
cluster of lattice elements j, k and i is introduced to construct G. Such a 
method is called the second approximation method in this paper. 

Similarly, by introducing other specification indices representing existences 
of polygonal and/or polyhedral correlations of lattice points and by adopting 
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probability variables corresponding to larger clusters of lattice elements in order 
to construct G, it becomes possible for this method to proceed to higher approxi- 
mations with various degrees. 

In the present paper, the probability variables YR 2j,k,i, etc., are called 
cluster variables, and the factor G in Eq. (A.1) is called combinatory factor 
Though this method has been called by various names, it seems most suitable 
to call cluster variation method, since in this method the term giving the 
maximum contribution to the partition function is selected from the summands 
of Eq. (A.1) by maximizing the contribution with respect to cluster variables, 


and the free energy of the assembly is estimated from this maximum term. 


B. Combinatory Factors for Higher Approximations 


When larger clusters of lattice elements are adopted so as to describe the 
structure of a lattice in greater detail, the method becomes tremendously 
difficult in the following points: 

a) Proper construction of the combinatory factor. 

b) Derivation of the free energy from the partition function. 

With regard to a), aside from the increasing difficulty in constructing the 
combinatory factor itself, it has been pointed out by several authors that the 
method could not necessarily be improved by merely adopting larger clusters. 
Proper construction does imply proper adoption of larger clusters as well. 

A general method for the proper construction of combinatory factors has 
been given by Kikuchi.» However, his method is somewhat difficult to 
manipulate. In the present paper, an alternative method is proposed. The 
method is quite systematic in constructing combinatory factors.. Furthermore, 
it not only makes the construction of combinatory factors easier, but also 
contributes much to clarifying the features of various degrees of approximation. 
This method is tentatively called pseudo-assembly method the details of which 
may be given in the next chapter. 


II. Pseudo-Assembly Method 


C. Construction of Combinatory Factor 
for One-dimensional Lattice 


In order to facilitate the understanding of the pseudo-assembly method®.” 
for constructing combinatory factors, and to explain the terminologies and 
abbreviated notations used in this method, an assembly of lattice elements on 
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a one-dimensional lattice is considered in the first place. 

When an assembly of N single-site lattice elements with » components on 
a one-dimensional lattice composed of N lattice points is considered, and the 
number of lattice elements of species j, j=1, 2,---, m, is written as Nj, it 
follows that 


SN = Ns (C.1) 


If the probability of finding a lattice element j on an arbitrarily selected 
lattice point is denoted by x;, and that of finding a pair cluster of lattice 
elements j and k on a pair of lattice points neighboring each other, i.e., on 
a lattice bond, by 4;,x, it follows that 


=NIN= > ys, (C.2) 
a 

Sere ie and > Yin = be (C.3) 
j=l jk=1 


In Fig. 1, a fragment of an assembly of lattice elements with two com- 
ponents is shown, and the explanation is given of the method for this assembly, 
but the formulation is made for a more general case of ~ components. It is 

—e—0—e—8—0—0—0—e— assumed that the assembly in (a) satisfies the 
iaytAssamblyt on ronecdi mensional lattice given sets of a; and y;,x. In such a case, a 
simpler expression “the assembly has a right 
distribution (r.d.) in y” is used. In (b), an 
$$ Oe @—0 O—O @ alternative representation of the same assembly 
(b) Assembly having rd in y is given so as to depict the assembly as a 
mode of configuration of pair clusters on lattice 
bonds. In this figure, every lattice element 
@—0 OO @ is represented doubly on the upper and lower 

io Gassabiy sides of each lattice point. 
The problem to construct the proper 
combinatory factor for a one-dimensional 

@®oe@eeo00e lattice is to find out the number of distinguish- 

"“@ 0o@eo0008 : 
の sea as able configurations G® of such an assembly 
of pair clusters as shown in (b). For this 
purpose the following procedure is devised. 


®oo0 @oo00e (1) The pair clusters are distributed on 
®@eeoeo0oe6eeoeé ; ; 
(e) P-assembly the lattice bonds independently to construct 


Fig. 1. an assembly of pair clusters having a right 
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distribution in y as shown in (c) which is called a Q-assembly. In the Q- 
assembly such a configuration as shown in (c) is necessarily included, and two 
different species of lattice elements occupy the same lattice point. The number 


of distinguishable configurations of the Q-assembly, denoted by QA is 
obviously given by 


Qn = NI/ I (ys.4N)! . (C.4) 
For the brevity of expression, the abbreviated notations similar to 
[point]: = L! (I 、 
[bond], = EY/ I (yseE)! , 


and relations similar to 


[point].z=((point]z) , 
[bond] qz=([bond]z)? , (C.5) 


will be used hereafter without remarks.* 
Adopting the abbreviated notation given above, Eq. (C.4) may be written 
as 


Q® =[bond]y . (C.6) 


(2) 2N lattice elements are distributed on the upper and lower sides of 
the lattice points to make an assembly of 2N lattice elements having right 
distribution only in % as shown in (e) which is called a P-assembly. The 
number of distinguishable configurations of the P-assembly is denoted by P™, 
and is given by 


P“=[point]z2y=([point]w)? . (Cay) 


(3) In the Q-assembly, there should be configurations in which two lattice 
elements located on any of the lattice points coincide with each other as 
shown in (b). The number of such configurations is surely G® in question. 

(4) Quite similarly, in the P-assembly, there should be configurations as 
shown in (d). The number of such configurations, denoted by G"), is apparently 
equal to the number of distinguishable configurations having a right distribution 
only in x. Hence, 


G™=[point]y . (C.8) 


* In the case where L is sufficiently large, it is possible to use Stirling’s formula, 
which gives Eq. (C.5). 
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(5) Writing the fraction of G® in Q® and that of GY in P™ as I’? and 


I) respectively, i.e., 
PM=G/Q), and PHV=GY/PH , (C.9) 
the comparison of the role of T'® with that of T'™ reveals an essential similarity 


between them: J" selects G configurations from Q configurations of the 
Q-assembly while TW selects GW configurations from PW configurations of the 


P-assembly. 
(6) Now it is assumed that the following identity holds: 
raal | (C.10) 
from which it follows that 
G%=GHS® | (C.11) 
with 
S%=Qe2/Pw | (CEI2) 


In the present case of a one-dimenesional lattice, G® is expressed by Eqs. 
(C.6) to (C.12) as 


G® =([point]y)~—'[bond]y , (C.13) 
or explicitly by Eq. (C.5) as 
G®=I(xjN)!/I(yi,4N)!. (C.14) 


In the present paper, the P-assembly is called the pseudo-assembly for the 
Q-assembly. Since the construction of the P-assembly plays an essential role 
in the method outlined above, it may be called the pseudo-assembly method. 

It can be shown in the following sections that the method is easily applied 
to two- or three-dimensional lattices. 

In the present case of a one-dimensional lattice, the exact expression of 
G for a circular lattice has been given by the present author* as 


G® =C,([point]v)~![bond]y , 
Cu = (MT yin 1s) , 
gj .k=1 j=l 
with 


Ti a= Lyin tt ek; 
ale inyaee (C.15) 


which holds for any N when Y;.x=Yk.s, i.e., the assembly is isotropic. 
By inspecting the terms in G™, the factor C, is found to be perfectly 


* See Ref. 6) where C, and Cs are given. C4 also conforms to Eq. (C.15), but the 
result is not published. For larger mn, practical calculations are too cumbersome to 
carry out, and the equation is the formal generalization of the forms for C2, Cs and Cy. 
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negligible in comparison with ([point]y)-![bond]y. Thus, it may be said that 
G” given in Eq. (C.13) by the pseudo-assembly method is exact from the 
standpoint of statistical thermodynamics,* and that the underlying assumption, 


given in Eq. (C.10), in the pseudo-assembly method is justified at least in the 
case of a one-dimensional lattice. 


D. Combinatory Factor for 1st Approximation 


Here the method developed in the previous section is applied to assemblies 
on a two- or three-dimensional lattice. As pointed out in Sec. A, hi=o is used 
as the sole specification index of lattice structure for the combinatory factor 
for the first approximation method. 

Consider an assembly of N single-site lattice elements on a lattice with NV 
lattice points and woN lattice bonds. In Fig. 2 (a), a small part of an assembly 


— 


(a) Original assembly having r.d. in y 


Bb: 58K 


(b) Assembly having r.d. in y {c) Q-assembly 
(d) Assembly having r.d. in z (e) P-assembly 
Rig. 


* The combinatory factor given by Eq. (C.13) for a one-dimensional lattice has been 
proposed by various authors, giving the same results as those obtained rigorously with 
the method of matrix. (See, for example, Ref. 1) of the present paper.) However, the 
direct comparison of G itself with the exact expression for any N and m does not 
seem to have been made. 
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with two components (7=2) is shown on a two-dimensional triangular lattice 
(w=3), and the explanation will be given for this case, but the formulation 
will be made for the general case of any m and w. In (b), an alternative 
representation of (a) is given, and in order to draw attention to the pair 
clusters, each of the lattice elements is written 2w-ply. (c), (d) and (e) correspond 
to those in Fig. 1, respectively. These figures themselves apparently show 
their constructions. 

The number of distinguishable configurations of the Q-assembly, Q, that of 
the P-assembly, PW, and that of the assembly in (d), GM, are given respectively by 


Q” =[bond].v=([bond]y)® , (D.1) 

PW =[point]:ux =({point]x)” , (D.2) 
and 

G“=[point]y . (D.3) 


By introducing these equations into Eqs. (C.11) and (C.12), the combinatory 
factor 
G.. =([point]»)'~2*([bond]) , (D.4) 


is obtained which has already been derived by Takagi‘) and others in rather 
complicated ways. For a one-dimensional lattice, by putting w=1, Eq. (D.4) 
is reduced to Eq. (C.13). 


EK. Combinatory Factor for 2nd Approximation 


In the first approximation method, since the combinatory factor is constructed 
by using as the sole specification index of lattice structure, there can be no 
distinction, for example, between a two-dimensional triangular lattice and a 
simple cubic lattice, both having the common coordination number, 20=6. In 
order to proceed to the second approximation method, it is at least necessary 
to indicate the fact that the former lattice has »N triangular cells by using » 
as the second specification index hs of the lattice structure. 

In Fig. 3(a), there is given a duplicate of Fig. 2, but in the present case, 
the assembly has a right distribution in triangular variables as well. A triangu- 
lar cluster variable is denoted by 2;,x,; which represents the probability of finding 
a triangular cluster of lattice elements j, k and 7 on a triangular cell. The 
combinatory factor for the second approximation method is written as G,® 
and v=2 in the present case of a two-dimensional triangular lattice, since the 
lattice has 2N triangular cells. Full explanation will be given for this case, 
but the formulation will be made in common with the case of a face-centered 
cubic lattice, i.e., »=6. An alternative representation of (a) is given in (b), in 
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(a) Original assembly having r.d. in z 


(d) Assembly having r.d. in y (e) P-assembly 


Fig. 3. 


which special attention is paid to triangular clusters by writing each of the 
pair clusters (3v/w)-ply. (c) is the Q-assembly for the present case. In (d), an 
assembly having a right distribution in y, but not in z, is given, but the pair 
clusters are written (3»/w)-play along the lattice bonds to provide convenience 
for imitating (b). (e) is the P-assembly. It must be kept in mind that in the 
P-assembly two lattice elements at the same corner of a triangular cell coincide 
with each other throughout the lattice. In other words, the pair clusters are 
properly connected at the corners of triangular cells. 
The number of distinguishable configurations of Q-assembly is given by 


Q =[triangle]yw=([triangle]y)” , (E.1) 
where the abbreviated notation 


[triangle];=L!/ Ii Gini)! , 
IVS + 


which is similar to Eq. (C.5), is used. 

Now, in the P-assembly, 3» pair clusters having a right distribution in y 
are properly connected at the corners of the triangular cells. Hence, the 
number of distinguishable configurations of the P-assembly, P, is obtained 
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by introducing @=1 and substituting 3»N for N in Eq. (D.4). 


P®=([point]syv)~'[bond]svy , 
=([point]w)-*”([bond]y)*” . (E.2) 


The number of distinguishable configurations of the assembly in (d) is 
apparently Go® itself. 

When written =G,/Q™, and assumed that the similar assumption as 
expressed by Eq. (C.10) is still valid for the present case, i.e., 


Pw, (E.3) 
it follows that 

Gy®=G.2S® , (E.4) 
with 

S8=Q3)/P 2 , (E.5) 


or, when Eq. (C.11) is introduced into Eq. (E.4), 
Gy®=GHUS@S@) 。 (E.6) 


The complete expression for Gy® is obtained, by combining Eqs. (D.4), 
(E.1), (E.2) and (E.5) with Eq. (E.4), as follows: 


Gy =([point])!-2°t**([bond]w)®-*"([triangle])” . (E.7) 


By putting @=3 and v=2, and w=6 and v=8, the combinatory factors of 
the second approximation method may be obtained for a two-dimensional 
triangular lattice and a face-centered cubic lattice, respectively. 

By the exactly similar procedure, the combinatory factors of the second 
approximation method for various kinds of lattices may be easily constructed. 
Many examples including those of the third approximation method have been 
given in Refs. 6) and 7). 


F. Selection Rule of Cluster and Degree of Approximation 


The generalization of Eqs. (C.11) and (E.6) may lead easily to the following 
equations: 


G=GHTTS®) ‘ 
{s} 
{s}=2, 3,° a (F.1) 


where G" represents the combinatory factor constructed in terms of r-lattice- 
point-ciuster, and 
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Ges lim G™ : (F.2) 
where G™ represents the rigorous combinatory factor with which the exact 
results would be obtained, in principle, by the cluster variation method. 

However, there still remains a problem how to choose the proper procession 
of s, i.e., the proper set {s}. 

In Bethe’s method, the cluster shown by (e) in Fig. 4 is used in the case 
of a two-dimensional square lattice. However, it has been shown that the 
adoption of Bethe’s cluster as well as those shown by (b) to (d) leads to the 


wo coo ob ob of. 


(a) (b) (c) (d) ie) 
(f) (g) | | 
(h) 
(i) 
Fig. 4. 


same results as obtained by using pair cluster shown by (a). Hence, the 
procession of s=2, 3, 4, 5 along the series of figures in the first row of Fig. 
4 is trivial in improving the degree of approximation. The present author 
and his collaborators®”) have found the similar situation that the adoption of 
clusters given as (g) and (h) does not improve the degree of approximation 
higher than the second approximation in which square cluster shown by (f) is 
adopted. Hence, the procession of s=4, 6, 8 along the series of figures in the 
second row is of no significance. 

On the other hand, it has been proved that the successive adoptions of 
clusters along the first column of Fig. 4, i.e., (a), (f), @, improve the degree 
of approximation step by step. Hence the procession of s should be taken 
along the first column, i.e., s=2 (pair), 4 (square), 9 (quadruple square), etc. 

Now it must be noted that the succession of clusters, (a), (f), (i),-+:, is 
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characterized by the following points: 

(1) Cluster (f) is closed with respect to cluster (a). 

(2) Cluster (i) is closed with respect to cluster (f). 

The similar situation has been found in other types of lattice, and the 
present author has generalized this situation 
and presented a rule for the selection of a) One-dimensional Lattice: OO 


the proper procession of s: b) Two-dimensional Lattice : 


: (i) Square lattic 
“Larger cluster should be closed with ilies oe: 
respect to smaller clusters.” 6—6 = Tal ; sei 7 


Some examples of proper succession of (i) Triangular lattice 


clusters are given in Fig. OE 
The rule given above seems to be some- oo -: LN y ee 


what empirical, but it has been justified more © Three-dimensional Lattice : 


rigorously by Kurata.® (i) Simple cubic lattice 
Hence, Eqs. (F.1) and (F.2) may be oo . r? : fee , ete, 


rewritten as 


(ii) Face-centered cubic lattice 


tA ae 


Gn =G II S$‘) ; 
a Fig. 5. 
{as}=a1, G2,°++, ar, (F.3) Proper Succession of Clusters 
and 
G@= lim G@” , (F.4) 


ro 


where a, represents the number of lattice elements in the s-th cluster in proper 
succession of clusters as shown in Fig. 5, for example. 

The cluster variation method using G'“” as the combinatory factor should 
be called the 7-th approximation method. Thus the degree of approximation 
which has been left quite ambiguous for the lack of criterion is assigned 
easily. 

The rate of improvement of the cluster variation method as 7 increases 
seems quite rapid, since the results of the second approximation method for a 
two-dimensional Ising lattice approach quite closely to the rigorous solution 
obtained by the matrix method. 

The present method of constructing combinatory factors has been studied 
by Hijmans and De Boer® with further examinations and applications. 


Bethe Lattice and Kramers-Wannier-Kikuchi Lattice 


The essential features of the first approximation and the second approxi- 
mation methods must be pointed out here, for they are very important in 


Some Contribution to Lattice Model Theory of Polymer Solutions 55 


treating polymer solutions by the present method. 

In the first approximation method, the lattice structure is specified only 
by its coordination number as pointed out frequently. Such an abstraction of 
the lattice structure is equivalent to replacing a real lattice with a hypothetical 
lattice having a dendritic structure shown in Fig. 6. The lattice may be 
called Bethe lattice, since Bethe’s approximation method, i.e., the first approxi- 
mation method, is exact for this lattice owing to the fact that any assembly 
on Bethe lattice does not contain any cluster 
closed with respect to pair clusters. This 
is the essential feature of the first approxi- 
mation method. 

Similarly, in the second approximation 
method, a real lattice is replaced with a 


lattice having only the smallest closed cells Bethe lattice with 
20 =4 


with respect to lattice bonds which are rahe 
ig. 6. 


connected with one another on every side 

of the closed cell as in the original lattice, nevertheless it has no larger closed 
cell of lattice points. For example, in the case of a two-dimensional triangular 
lattice, the lattice is replaced with a hypothetical lattice having only triangular 
cells connected side by side, nevertheless it lacks in hexagonal correlation of 
lattice points. Such a lattice may be called Kramers- Wannier-Kikuchi lattice, 
since it has been found by Kikuchi® that Kramers and Wannier’s approximate 


method is equivalent to the second approximation method. 


G. Combinatory Factor for Intermediate Approximation 


It may be interesting to consider the construction of the combinatory 
factor for the cluster variation method with intermediate degree of approxima- 
tion.* The combinatory factor for the (1+g)-th approximation method, (0<q<1), 
for triangular lattices may be constructed by applying the pseudo-assembly 


method only to gyN triangular cells. The combinatory factor, G\), may be 


given as 


GA =([point]w)*?°t* »([bond]w)®~*4"([triangle]w)2” , (G.1) 


which may be reduced to Eqs. (D.4) and (E.7) by putting g=0 and q=1, 


respectively. 
When gq satisfies specially the condition: 


* Somewhat generalized treatments of combinatory factors for intermediate approxi- 
mations are given in Ref. 7). 
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qg=o/(3y) , (G.2) 
the combinatory factor takes a very simple form: 
G),=([point]w)'~°((triangle]w)°/* (G.3) 


This is essentially equivalent to the combinatory factor given by Yang, 
Li and Hill.'® 


(a) Original lattice (b) Yang-Li-Hill lattice 


Fig. 7. 


The lattice on which Eq. (G.3) holds rigorously is such a lattice as shown 
in Fig. 7 which may be called Yang-Li-Hill lattice. This lattice has no closed 
cell larger than a triangular cell, and triangular cells are connected with one 
another at their corners so as to satisfy the coordination number of the original 
lattice. 


Ill. Chain Polymer Solutions 


Since the earlier development of the lattice model theory of polymer solu- 
tions which owes much to pioneering works of Flory and Huggins, a large 
number of improved theories have been proposed by various authors.!” These 
theories deal chiefly with the combinatory factor or geometrical effect of chain 
polymers. However, there are other two effects, energetical effect and structural 
effect. The energetical effect is due to the difference in interaction energies 
between two neighboring lattice elements. This kind of effect is exactly the 
same asin monomer solutions. On the other hand, as pointed out by Prigogine, '”) 
the structural effect is specific to polymer solutions and due to the difference 
in external degrees of freedom of various lattice elements. However, one of 
the main defects of existing theories seems to rise from the incomplete estima- 
tion of the geometrical effect of chain polymer molecules. 

It may be seen easily that the majority of existing theories correspond to 
the first approximation method which uses a Bethe lattice shown in Fig. 6 of 
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Sec. F, since in these theories the lattice is specified only by its coordination 
number. In order to make the feature of existing theories clearer, consider, 
for example, an assembly of chain polymer molecules on a two-dimensional 
triangular lattice. Three consecutive segments, labeled by j, j++1 and 7+2,/ of 
a certain polymer chain may occupy three consecutive lattice points in various 
ways. When the three segments occupy a triangular cell of the lattice, the 
segments j and 7+2 may be found respectively on two lattice points neighbor- 
ing each other. Such an intramolecular correlation of the segments j and j+2 
is called an inner correlation of a polymer molecule. In existing theories, the 
triangular lattice is replaced implicitly by a corresponding Bethe lattice which 
is lacking in triangular cells. Hence, such an inner correlation cannot be 
transcribed onto the Bethe lattice. Furthermore, on the Bethe lattice, once 
two separate chain polymer molecules approach each other, and two segments, 
each belonging to a separate polymer molecule, occupy two neighboring lattice 
points, the other parts of the two polymer molecules can never approach each 
other to occupy another pair of neighboring lattice points as seen clearly 
from Fig. 6. In brief, intermolecular pair correlation of two separate chain 
polymer molecules can be taken into account once and only once. Such 
situations as neglect all of the inner correlations and as limit the intermolecular 
correlations of two separate polymer molecules to only once are surely drastic 
oversimplification of the real situations on the original triangular lattice except 
one case where the polymer molecules are stretched rigidly or rod-like. 

In order to take the shape of a polymer chain into account, it may be 
necessary to proceed to higher approximation methods in which closed correla- 
tions of lattice points are taken into account. 

In the following sections of the present paper, an improved lattice model 
theory of chain polymer solutions is given in order to calculate the second 
osmotic virial coefficient in terms of shape factors which define the shape of a 
polymer chain. Though the method of the calculation may be applied quite 
generally to assemblies on various kinds of lattices, the practical] calculation 
is carried out in the case of triangular lattices and the greater part of numerical 
results for the second osmotic virial coefficient are given in a more limited case 
of a two-dimensional triangular lattice for the purpose of exemplifying this 
general method. 

Since the original type of this improved theory has been given in Ref. 13), 
it may be convenient for some readers to read the following sections in close 
conjunction with Parts I and II of the referred paper. 
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H. Basic Formulation 


Consider a lattice model polymer solution composed of No single-site solvent 
molecules labeled by j=0 and N, p-site chain polymer molecules. Each of the 
polymer molecules consists of p segments which are connected by p—1 chemical 
bonds and are labeled by j=1, 2,---f from the head to the tail. 

By writing the total number of solvent molecules and segments, which is 
equal to the number of lattice points, by N, and the volume fraction of polymer 
in the solution by 1, it follows that 

N=N.+fN> , (H.1) 
and 
v=pN,/N . (H.2) 

In Fig. 8, a small part of the lattice model 
polymer solution is given on a two-dimensional 
triangular lattice, and here the _ solvent 
molecules are omitted in order to simplify the 
figure. However, the formulation is made as 
far as possible for a triangular lattice with 
oN lattice bonds and vN triangular cells. A 
two-dimensional triangular lattice corresponds 
to w=3 and v=2, while a face-centered cubic 


lattice corresponds to w=6 and »=8. 


Fig. 8. 


The chemical bond connecting a segment 
j with a segment j+1 is shown by a thick solid line with an arrow indicating 
the sense of the chemical bond as shown in the figure. A pair of consecutive 
segments j and j+1 connected by a chemical bond is written as j>,<j+1 
symbolically. The inner correlation of the segments j and j+2 is shown by 
a thin broken line with an arrow indicating the sense of this correlation as 
shown in the figure. The broken line is called an inner correlation bond. A 
pair of segments j and j+2 connected with an inner correlation bond is writ- 
ten as j>,Cj+2 symbolically. 


Restriction on Intramolecular Correlations 


In Fig. 9, is given another part of the same polymer 
= >< 


solution in which various kinds of intramolecular 

correlations are shown. However, the following restric- NN 
tion is imposed on the occurrence of intramolecular Ny 
correlations: lex 


“Intramolecular correlations between two segments j 
and j+s with s=3 are ignored, i.e., prohibited.” Fig. 9. 
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Hence, the existing intramolecular correlation is only the correlation between 
segments j and j+2 which has been called an inner correlation as mentioned 
previously. 

The restriction given above is partly intrinsic in the present theory and 
partly for the convenience of simplifying the formulation. However, the detailed 
consideration of this restriction will be given later in Sec. M. 


Shape Factors 


In order to specify the shape of a polymer chain on the original lattice, 
shape factors 6B;®, s=1, 2,---, and j=1, 2,---, p, are introduced for each of 
the segments of a polymer molecule in the following ways: 

6;: Probability of being bent z/3 at the segment j of three consective 
segments j—1, 7 and j+1. Hence, the probability of being not bent 
z/3 is 1—B;. 

6;: Probability of being bent 27/3 at the segment 7 of the three consecu- 
tive segments not bent z/3 at the segment 7. Hence, the probability 
of being bent 27/3 at the segment j of the three consecutive segments 
is (1—f;'Y)8;”, and that of being kept stretched at the segment 7 is 
(1— pj) — 8). 

8;: Probability of finding a bending in cis form of four consecutive 
segments j—2, j—1, 7 and j+1 after having been bent 27/3 twice at 
the segments j—1 and 7. Hence, the probability of being bent in cis 
form of the four consecutive segments is 


(1-6 DBP LBL, 
and that in trans form is 


(1p )8%, BC . 


For the full description of the shape of a polymer chain on the original 
lattice, a large number of the other shape factors are necessary. However, 
the three kinds of shape factors given above are sufficient for the present 
formulation. 

When 


B;=0 or 1, for all s and 7, 
the polymer chain is rigid and takes a definite shape. On the contrary, when 
0<f;"<1, at least for some s and J, 


the polymer chain is flexible. Some examples of the relationship between 
shape factors and the shape of a polymer chain are given in Tables 5 and 6 


of Sec. L. 


60 T. WATARI 


Cluster Variables 


Consider an arbitrarily selected triangular cell abc composed of three lattice 
points a, 6 and c on which various kinds of cluster variables are introduced 
as follows: 

xj: Probability of finding a lattice element j on a lattice point, say a, of 
the triangular cell abc. ao represents the said probability of a solvent 
molecule, while &;, j=1, 2,---, p, represents that of a segment /. 

Yj>,<j+1: Probability of finding a pair of segments 7 and j+1 connected by 
a chemical bond on a lattice bond, say ab, of the triangular cell abc. 
Such a pair of segments has been written j7>,<j+1 symbolically. 

Yi>.cj+2: Probability of finding a pair of segments j and j+2 connected by 
an inner correlation bond on a lattice bond, say ab, of the triangular 
cell abc. Such a pair of segments has been written j>,Cj+2 sym- 
bolically. 

Z>j>,<j+1>,<j4+2¢: Probability of finding a set of three consecutive segments 
j, j+1 and j+2 on the triangular cell abc, j on a, j7+1 on b and j+2 
on c. Hence, the three consecutive segments are bent z/3 at the seg- 
ment j+1. 

These variables are connected with the volume fraction, v, of polymer by 
the following equation: 


aj=(lp), j=l, 2,---, Pp, 


=l-v, a 
Yi>.<i+1=(1/2e)(v/p) , gly eee Pl, 
=0 j=0, p, 
Y j>,c5+2= (Bj t1/2a)(v/p) , be 前 な ここ 
= j=0, pI, 2, 
255>,<541>,<542c =(By 41/62) (v/p) ; t food Wie’ as Oe Be ea a 
=0, j=0, p—1, p. (H.3) 


Eq. (H.3) may be obtained by the following considerations: The probability 
of finding the segment j on the first lattice point a of the triangular cell abc 
is apparently x;=(v/p), and that of finding the segment j+1 connected with 
the segment j by a chemical bond on the second lattice point 6 is (1/2w). 
Hence the probability of finding 7>,<j+1 on the lattice bond ab is (v/p)-(1/2@) 
which is the second equation. Consider three consecutive segments j, j+1 and 
j+2 of which j>,<j+1 has occupied the lattice bond ab. The probability of 
being bent z/3 on the lattice point b at the segment j+1 of the three consecutive 
segments is 6A,, and that of finding the segment j+2 on the lattic point c is 
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(1/4), where 2 is the number of triangular cells which are connected with one 
another at a common lattice bond, say ab, and is given by 2=(3»/w). Hence, 
the total probability of finding the three consecutive segments on the triangular 
cell abc, j7 on a, j+1 on b and j+2 onc, is Yj>,<441° BS 1:(@/3v) which gives 
the fourth equation. Since the pair of segments j つ , こ j+2 found on the lattice 
bond ac of the triangular cell abc appears 2 times on 2 different triangular 
cells, the triangular cell abc inclusive, the probability of finding 75,cj+2 on 
a bond of an arbitrarily selected triangular cell is 25;s.<j415,<j420:4 which 
leads to the third equation. 
The following cluster variables are also used: 

Yj,x: Probability of finding a pair cluster of lattice elements j and k on a 
lattice bond, say ab, of the triangular cell abc, 7 on a and k on b, each 
of the lattice elements belonging to a separate molecule. 

2j,k,i. Probability of finding a triangular cluster of lattice elements 7, k and 
z on the triangular cell abc, j on a, k on 6b and i on c, each of the 
lattice elements belonging to a separate molecule. 

Zj>,<j41,k: Probability of finding a triangular cluster of a pair of segments 
j>,<j+1 and a lattice element k on the triangular cell abc, j7>,<j+1 
on ab and k onc, each of j>,<j+1 and k belonging to a separate 
molecule. 

Zj>,cj+2,k: Probability of finding Table 1. 

a triangular cluster of a pair 
of segments j>,Cj+2 and 
a lattice element k on the 
triangular cell abc, j>,Cj+2 
on ab and k on c, each of 
j>,Cj+2 and & belonging 


to a separate molecule. 
The types of cluster shown by the subindices of cluster variable are given 
in Table 1 which may help toward the understanding of these cluster variables. 


Second Subindexing of Cluster Variables 


When special attention is paid to a segment j locating at a lattice point, 
say a, of a triangular cell, say abc, of the original lattice, there will be found 
two branches j, j—1, j—2,---, 1 and j, j+1, j+2,---, p of a polymer chain 
starting the segment j. The configuration of these two branch chains on the 
lattice may be found in various modes with respect to the triangular cell atc. 
Each of the modes of configurations is specified by an index s. When the 
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configuration of the two branch chains is found in the s-th mode, the segment 
j is written as j;s symbolically, The probobility of finding a segment Fes 
on a lattice point of a triangular cell is written as x;,, by putting the index 
s to 2; as the second subindex. Similarly, y;,x;s,. represents the probability 
of finding a pair cluster of two independent segments j ; s and k ; ¢ on a lattice 
bond of a triangular cell. When k=0, for example, k represents solvent 
molecule, and the second subindex ¢ for k is unnecessary. However, in order 
to keep the symmetry of the second subindexing, y;.0;s,.0 is used in this case. 
Zj,k.i;rs,tu,ow Yepresents the probability of finding a triangular cluster of three 
independent segments j, k and i composed of j,k; s,¢ on ah, k,i;u,v on be 
and i,j;w,r on ca. In this case, the same segment j, for example, is 
subindexed doubly by 7 and s for the purpose of making the decomposition of 
j,k, i; rs, tu, mw into i,k; s,t, k,i;u,v and i,j; w,7r easier. Similar second 
subindexings are made on the other cluster variables. 

Though such heavy subindexings seem to be quite cumbersome for the 
formulation of the present theory, they make it possible for the theory to give 
a rigorous formulation as will be seen later. 

Since the shape of a polymer chain is fixed by shape factors, {f;‘}, 
definitely or in average, the modes of configurations of two branch chains occur 
in definite fractions. Thus, «;.;, for example, is a definite fraction of xz; which 
has been connected with the volume fraction v of polymer by Eq. (H.3). It 
follows, therefore, that 


Li 32=Pj:0 , (H.4) 


where ¢;;;, is a factor which will be determined by giving {6;®} for a given 
lattice. Quite similarly, 


Yi>,<541; re=Pj>,<j413 red , 
Y5>,6542:r.0=Pjr,cjt+2; 7.80 , 


Z>j>,<j+i>,<j+2c 7 7rs,tu,vw =$355>,<j41>,<j+2c srs,tu,vwV , (H.5) 


where ¢’s are factors similar to ¢;,.. 


Partition Function and Free Energy 


The following formulation is given for the (1+q)-th approximation method. 
When q=0, q=(@/3v) and g=1, the formulation corresponds to the first (Bethe- 
Fowler-Takagi), the [1+(/3v)]-th (Yang-Li-Hill) and the second (Kramers- 
Wannier-Kikuchi) approximation methods, respectively, as mentioned in Sec. D. 

A term in the summands of Eq. (A.1) for the configurational partitior 
function is written as 
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Z=Gyy exp [—E/kT] , (H.6) 


(3) : : 
where G;, is the combinatory factor which has been formally given by Eq. 


(G.1) with abbreviated notations. 


The abbreviated notations are given, in the present case, as follows: 


[point]y =WN!/II(a;.,N)! , 
[bond] 7 =N L/II(y3.x : rsV) ! [Yi> Si ; rsV)! (Yjr,c542 ; re)! ]? ’ 
[triangle]ly=N MZ : pein aN) ! [(2j>.<j41.6 ;7rs,tu,vw N)! r 


“4 (275,67 42,% ‘ panei On Sauk pee eae ‘ V0) Ny ’ (H.7) 


where the products are made over all of the existing cluster variables of types 
written after the product signs. The powers 2 and 6 raised to the terms be- 
tween square brackets come from the fact that the clusters of lattice elements 
corresponding to the cluster variables appearing in these bracketed terms appear 
on a lattice bond and on a triangular cell in 2 and 6 ways, respectively. 

The interaction energy, written as E in Eq. (H.6), is given as follows: 


p 
E=oN(> 565,59 + Y;,4Ae;,n) , 
rE 


p 
j,k=0 
Ae; n=; ,h—(€7,5 +€x,%)/2 , (H.8) 
where ¢;,z represents the interaction energy contribution of a pair cluster of 
independent lattice elements 7 and k. 


Instead of maximizing Z with respect to independent cluster variables, g 
given below is minimized. 


u=—(RT/N)1nZ 
=—(kT/N)1nG) +E. (H.9) 


The minimizations give a group of equations for the so-called quasi-chemical 
equilibria between clusters of lattice elements which connect cluster variables 
one another. The procedure of the minimizations is given later together with 
the equations for quasi-chemical equilibria. 

With the equations for quasi-chemical equilibria, and Stirling’s formula, Eq. 
(H.9) takes the following simple form after somewhat cumbersome rewritings: 


L=kT > [1 —204+3q»)x;; nv; ;s 
+(w—3qv)Y5,5 tne 1ny;,3 37,8 


+QZj,5,3; retanaer LOS sees ; tee! 


at > WX5E5,5 , (H.10) 
f= 


“where the first summation is made over all of the cluster variables written 


after the summation sign. 


64 T. WATARI 


# given above is now regarded as the free energy per lattice point, and 
the chemical potential per molecule of the solvent, go, is calculated from Eq. 
(H.10) by 


po=ONp/ONo=p—v(Op/dv) . 
The calculation gives 


(p40 — po) /RT =(1—2@+3qv)1na0+(@—3qv)1nYyo,o 
+qvl1nZo,0,0 , (H.11) 


where the abbreviation, 
fo° =E0,0 = 


is used. 


Equations of Quasi-Chemical Equilibria 


Upon minimizing g given by Eq. (H.9) with respect to cluster variables, 
the choice of independent variables is made as follows: 


25,k,i; r8,tu,vw with the exception of 
27,0,0; rs,00,00 » 

2j>,<g+1,.k ; re.tu.ow with the exception of 
Zy>,<j+1,0; rs,tu.00 5 

2j>,CH+2,k 3 re,tu,vw with the exception of 


Z2j>,Cj+2,0; rs,tu,00 » 
and 


v. 


For the minimization with respect to 2;,x,0; rs,tu.00, Which is a special case 
Of 2;,4,i;re,tu.ow With i=0, the following relations are used: (2) -ale re d= 


92j,0,0 ; r2,00,00/02=—1 , 

20, k,0 ; 00,tu,00/0Z= —1 ’ 

020,0,0/0z=1 , 

Yj.k ; 0,t/0Z=1 , 

9Y5,0 ; »,0/0z=—1, 

0Yo,k ;0,1/02=—1 , 

Yo,o /Az=1 , (H.12) 


which may be obtained from the following equations connecting these variables 
with Zj,k,0; rs,tu,00 : 


Uj ses =2j,0,0; r8,00,00+Zj,k,0; rg,tu.00 十 … 


a 
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Lk st =Zo,k,0 ; 00,tu,00+2Zj,4,0; rs,tu,oote*** 


vo =Z0,0,0 +20,k,0 ; 00,tu,00 十 … , 


Yj,k; 5,4=25.k,0 ; ratuooters: , 


Yj,0: s,0=25,0,0 ; rs,00,00 十 ・・ 
Yo,k ; 0,t=20,k.9 ; 00.tu,00 二 **** , 
Xo =Yo.0o+Yo.k:o2et°°°* , 
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(H.13) 


and Eqs. (H.3) and (H.4) showing x’s are definite fractions of volume fraction 


v of polymer for a given set of shape factors. 


By remembering Eq. (H.12), the minimization gives the following equation 


of quasi-chemical equilibrium: 


FE ; 0.tYo.o)/(Ys.0; 2,0Yo,k ; 0,1? °%: 
-[(23,%.0 ; rs,tw,.00Z0,0,0)/(2;,0,0 ; 18,00,0020,k,0 ; 00,tu,00) 


SWE leh Eble 


1820 


Putting g=0, 
(Y5,% : s,tYo,0)/(Ys,0; s,0Yo,r 0 = Hd HHA; 
hence, 


(Z;,%,0 ; rs,tu,00Z0,0,0)/(2;,0,0 ; rs,00,00 200,k,0 ; 00, tw, 00) 
=[(Aj,0Ho,%/Hj,x)'” , 
where, 
H;.x=exp [(24e;,x)/(RT)] , 
and, by Eq. (H.8), 
guia. 
By exactly similar procedures, the minimizations with respect to 
Zj>,<j+1,k 3 7rs,tu,vw 


and 


2j>,Ci+2,k ; rs.tu,vw 


give the following equations, respectively, corresponding to Eq. (H.16): 


(Zj>,<p410k ; ra,tu,vw20,0,0)/(Zj>.<541.k ; rs,tw,0020,0,k ; 00,00, vw) 


=[Ho,3Hj41,0H oH, «/(He,iHi+1,%)]'” > 
(2)35,c5+2,% ; rs,tu,vw20,0,0)/(Zj>,c4+2,0 ; re,tw,0020,0,k ; 00,00, vw) 


= [Ho,jHs+2.0Lle,oHo,«/(He,iHi+2,%)|'”? . 


Osmotic Virial Coefficients 


The osmotic pressure z is given, in the present case, by 


(H.14) 


(H.15) 


(H.16) 


(H.17) 


(H.18) 


(H.19) 


(H.20) 
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r= —(po— fo°)/Vo ’ 


where w is the molecular volume of the solvent of a polymer solution. With 
Eq. (H.11), the above equation becomes 


(xv0/kT) = —(1—2w+3qv)1nao—(w—3qv)1nyo,o 
—qv1nZo,0,0 , (H.21) 


which may be expanded in the power series of the volume fraction v of polymer 


as follows: 
(xv0/RT)= >: dav” , (H.22) 
n=1 
in which bs is the m-th osmotic virial coefficient. 
Writing quite formally, 


Yo,o0 =1+a0,0v+bo,0v? +O(v?) , 
Z0,0,o=1+Ao,0,070+Bo,0,0v? +O(v*) , (H-23) 


and using Eq. (H.3) for xo, the firist and the second virial coefficients may be 
obtained from Eq. (H.22) as follows: 


bi=(1—20@+3qv)—(w—3qv)ao,o—qvAo.o,0 , (H.24) 
b2=(1/2)(1—20+3qv) —(@—3qv) [bo ,o—(1/2)a5,0] 
—q»[Bo,o,o—(1/2)Ad,0,0] . (H.25) 


Non-Athermal and Athermal Solutions 


In the following formulation, the segments of polymer molecules are assumed 
to possess the same chemical and physical natures. It follows, therefore, that 


Ae;, n=Ae , j>0 and k=0, or j=0 and k>0, 
=O, 4, k=0..0r 7, k>0; (H.26) 
and 
Ay.= HH, j>0 and k=0, or j=0 and k>0, 
==] j,kR=0,0r- 3 k>0 . (H.27) 
When 
4e+0, and H+1, (Non-athermal) (H.28) 


the polymer solution is called non-athermal. 
When 


46=0 , and A=1 5, (Athermal) (H.29) 


the polymer solution is called athermal. 
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I. First (Bethe-Fowler-Takagi) Approximation 


The first, i.e., Bethe-Fowler-Takagi, approximation method corresponds to 
the special case g=0 of the preceding formulation. Since the first approxima- 
tion method for polymer solutions has been studied extensively by many authors, 
none of new results may be obtained from the present formulation. However, 
in order to render help to the understanding of the method of calculating 
osmotic virial coefficients by higher approximation methods, and for the con- 
venience of comparing results obtained by various approximation methods with 
one another, the calculation is began with the use of the first approximation 
method. 

In the first approximation method, the original lattice is replaced by a 
Bethe lattice shown in Fig. 6. Since Bethe lattice has no triangular cells, the 
basic formulation for this approximation method is obtained by putting formally 
£;"=0 as well as q=0. 

The equation of quasi-chemical equilibrium for the present case is given 
by Eq. (H.15) which may be written with Eq. (H.27) as follows: 


CE: ; s,tYo,0)/(Yi,0 ;8,0Yo,k ; 4) =H, : (1.1) 
Since, 
co 
=Yo,0 = >> > Yo, sj $0,8 3 
j>0 s 
Pjis =Y5.0; 8.0 + Sy Daag es ks 
k>0 t 
it follows that 
=== Dy ee Yj,k i; 8.t > 
j>0 g,k>t st 
=1+0O(v) , 
Y5,0;8,0—Vj;8— Oa Sah Yj,k; 3,t » 
k>0 t 
=Zjia +O), (1.2) 
where the summation sign 1 represents that the summation must be 


extended over all of segments j; s which may form pair clusters 7,0; s, 0.. 
By introducing these equations into Eq. (1.1), the following equation may 


be obtained: 
Yj 5 9,t= HU; ; he; :+O(v') . 


With this result, Eq. (1.2) for yo,o becomes 
yoo=l—v— DF) a5 yet H Dy gg ate FOO) , (1.3) 
3>0 ae 
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where the summation sign [34 represents that the summation must be made 
over all of segments j;s and k;¢ which may form pair clusters j,k; s,t. 
In order to carry out these summations, a new cluster variable &; is in- 
troduced as follows: 
€;: Probability of finding a lattice element j on a lattice point, say a, at an 
end of a lattice bond, say ab, with the lattice bond ab unoccupied by a 
polymer chain. 
The new cluster variable is connected with variables % and y by the fol- 


lowing relations: 


j= 09 — Yas <5 ta bei es 


which becomes with Eq. (H.3) 


€;=1—0 5 7=0 ; 
=[1—(1/2@)](v/p) , j=1,?P, 
=[1—(1/o)](v/p) , j=2, Dy" ty p=1 , 
and 
3 y=1—(aloyv , 
with 


a=1—(1/p) . (1.4) 
Using this new variable, the summations appearing in Eq. (1.3) may be 
expressed as 


>; 19, Vg 5 — 2 t= 24 to 
j I= 


j>0 


SU lar; gk es=( 2 £5) —2bo © £36) 
j= j= 


j,k>0 
and Eq. (1.3) becomes 


Yo0=1+[(a/w)—2]v+A [(a/w)—1]?v? +O(v') , 
from which 
a,0=(a/w)—2 , 
bo,o=H [(a/w)—1]? . 


When these ao,0 and bh.o are introduced into Eqs. (H.24) and (H.25), with 
q=0 in mind, it finally follows that 


bi=1/p ? 
b2=(1/2)—o[1—(1/2){(a/w)—2}* +H {(a/o)—1}] . (1.5) 


When athermal, i.e., H=1, Eq. (1.5) may be reduced to 
b:=(1/2)—(1/2)(a?/o) . (Athermal) (I.6) 
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Though this is the well-known result for the second osmotic virial coefficient 
of athermal polymer solutions obtained by various authors, it is derived here 
as a special case of Eq. (I.5) for non-athermal solutions. 

The calculation of varial coefficient has been carried out as the special 
case q=0 of the (1+q)-th approximation method for a triangular lattice, It 
must be noted, however, that the final result of virial coefficient is quite general 
.in any Bethe lattice, since a lattice has lost its original structure during the 
transcription onto a Bethe lattice. 


J. Further Formulation for (1+q)-th Approximation 


The formulation given here runs parallel to Eqs. (I.1) to (1.3) of the preced- 
ing section starting equations of quasi-chemical equilibria given by Eqs. (H.16), 
(H.19) and (H.20) for 0<gS1, which may be written with Eq. (H.27) as follows: 


(Z3,k,0; rs,tu,00Z0,0,0)/(2Zj,0,0 ; r#,00,0020,k,0 ; 00,tu,00 )=H , (J.1) 
(25>, <5-+1,k ; rs,tu,vwZ0,0,0)/(Zj>,<5+1,0 ; r8,tu,00Z0,0,k ; 00,00,0w) =H? , (J.2) 
(Zj3,c5+2,k ; r8,tu,vw%o,0,0)/(Zj>,c4+2,0 ; r8,tu,00Z0,0,k ; 00,00,0w) =H? . (J.3) 
Since, 
to—1—v 


=Z0,0,0+ >) Dd! Zo,k.0; 00,¢u,00-+ 24 Ds 2v,0,¢ ; 00,00, 00 
k>0 t i>0 v 

+2 >) Dy (Ze) <e4210; r8,20,00 tb 285,Ck+9;0 : 40,0000) » (J.4) 

Ek s,t 
Pj ; e=2j,0,0; rs,00,00+ 2s 24 2j,k,0; re,tu,o0b 3 2 の 0 
k>0 t i>0 v 

Ee Dy D3, (Zk>,<k41,5 ; tu,vw,rs 十 2Z4 っ , こ k+2,5 : Bg eae ’ (J.5) 

ku,v 


it follows that 
Zo,0,0=1+O(v) 5 
2j,0,0 ; rs,00,00 — Wj ; stO(v?) . (J.6) 
By introducing these equation into Eq. (J.1), the following equation may 


be obtained: 


2j,k,0; ra,tu,00— 145 ; 3k ; -+O(v') . NY 
Furthermore, since 


Yj>,<j+li;s,t —2j>,<G+1,0; rs,tu,00 


+ Dy, Dy Zgj>,<g+1.k ; rs,tu,vw » (J.8) 
k>0 v 


it follows that 
Zj>,<j41,0; r8,tu,00— Yi>,.<j+1; 8.t +O(v?) : (J.9) 
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When this equation and Eq. (J.6) are introduced into Eq. (J.2), the following 
equation may be obtained: 
る jj>,< く 3 十 1k: に 8,t Lk: »+O(v') ’ (J.10) 
Quite similarly, from 


Yj>,cj4+2; s,t— 2jd,Cj+2,0; rs,tu.00 


中 DY 3 Zj>,Cj+2.k ; rs,tu,ow (J.11) 
k>0 v 
and Ea. (J.3), 
2j>,.C5j+2,k 3 rs,tu,vw = H*Yj3,¢5+2 5 8,t Uk ; »tO(v’) ’ (22) 


By combining Eq. (J.5) to (J.12) with Eq. (J.4), Eq. (J.4) for zo,0.0 becomes 
20,0,.0=1—(v+2 3 [3.0,0] Vj;s 
j>0 
doy By [3> < 0 jis,t 
j>0 
+2 BUI: peared 
j>0 
+3H( 3 Bhp;, ess 
jk>0 
+2H 3 > TL; <j Hl; atCk 9 
3, k>0 
+2H = C4>,CI+2 Ky — 549; sk v) +O(V®) 4 (J.13) 
j,k>0 


where specially devised summation signs similar to those used in Eq. (1.3) are 
used, meanings of which will be given in detail in the paragraph following 


Eq. (J.19). 
Since, 


Yo,0o=Zo0,0,0+ p> Z0,0,k ;00,00,0w » (J.14) 


it follows, when Eqs. (J.5) to (J.13) are combined with Eq. (J.14), that 
Yo,o=l—(wt+ Dy Fd Oe Sg 
j>o0 
+2 > る Re 
j>0 
a at LI>+SI4301 ee 
j>0 
+H ( 2 C5. kOe. est 
g,k>0 
+4H DS) > SH ys, ju tbhio 
j, k>0 


+4H > Cj>,C9+2, Ky j5,cj+9:0t€k 0) +O(v') ° (J.15) 


a, k>0 
From Eqs. (J.15) and (J.13), 


ano j,0,0]p. j>,<j+1, . 
w= + 2 Wi sat 2 Bd? SEE seg cat 
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Poyeo2. olan. ? 
+2 2 ¥ YjB,;ci +t;s, 0 ’ (J.16) 


bo=H( SS CG, Ole, vest 


j,k>0 


ES stl ily ey ca) tees 


i, k>0 
j>,CI+2,k n 
4 es 0 (J.17) 
Ao,o,o=—(v+2 3 [5,0,0]gp He > CI> <I My oS csaa sat 
j>0 j>0 
Bet date lamest ie, (J.18) 
j>0 


Boo 0o=3AH( > (Fg, . are. ¢ 
j,k>0 


[j>.<j+1, k 
+2H SY EO Y ciara tls « 
IK 


+2H 2S CI5,C942, Hay K-51 0.620: 0)/V® . (J.19) 


3,k>0 


In Eqs. (J.13) to (J.19), specially devised summation signs are used in such 
ways that the sign 7.0.0] represents a summation over all of segments j;s 
which may form triangular clusters 7,0,0;7s,00,00, the sign S{%>-<#+10] a 
summation over all of clusters j7>,<j+1;s,f which may form triangular 
clusters j>,<j+1,0;7rs,tu,00, and the sign S{'/>-<i+2°1 a summation over 
all of clusters j>, Cj+2;5s,t which may form triangular clusters j7>, Cj+2, 
0; vs, tu, 00. Furthermore, the sign £0 represents a summation over all of 
pairs of segments 7; s and k; t which may form triangular clusters j,k, 0; 7s, tu, 
00, the sign >{'’>:<7+1-4] a summation over all of pairs each of which consisting 
of a cluster j>,<j+1;s,t and a segment k; v which may form a triangular 
cluster j>,<j+1,k; 7s, tu, vw, and the sign >.<3+%k] a summation over all 
of pairs each of which consisting of a cluster j>,Cj+2;s,¢ and a segment 
k;v which may form a triangular cluster j>,Cj+2,k; 71s, tu, vw. 

In order to carry out these summations, a group of new cluster variables 
are introduced as follows: 

€;: Probability of finding a lattice element j on a lattice point, say a, ofa 
triangular cell, say abc, with the lattice bonds ab and ac unoccupied by 
polymer chains and/or an inner correlation bond. 

1;>.<j+1 : Probability of finding a cluster j>,<j+1 on a lattice bond, say 
ab, of a triangular cell, say abc, with the lattice bond ac and bc unoc- 
cupied by a polymer chain and an Table 2. 
inner correlation bond. 

Nj>,cj+2 : Probability of finding a cluster 
j>,cj+2 on a lattice bond, say ab, 
of a triangular cell, say abc, with the 
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lattice bonds ac and bc unoccupied by a polymer chain. 
These new variables are shown diagrammatically in Table 2 in order to 
make them easier to understand. 
The relations corresponding to Eq. (1.4) for these variables will be given 
in the following sections. 


K. (1+/3v)-th (Yang-Li-Hill) Approximation 


The (1+/3v)-th, i.e., Yang-Li-Hill, approximation method corresponds to 
the special case g=o/(3v) of the preceding formulation. In this approximation 
method, the original lattice is replaced by a Yang-Li-Hill lattice shown in Fig. 
7. Once the configuration of polymer chains has been transcribed from the 
original lattice onto the corresponding Yang-Li-Hill lattice, it may easily be 
seen that, during this transcription, all of clusters j5,Ccj+2,0 and jD,Ccj 
+2,k have disappeared either by the disappearance of triangular cells on 
which these clusters were found or by the breakdown of inner correlation 
bonds involved in these clusters. Hence the basic formulation for this 
approximation method is obtained by putting formally 2j5,<¢j42,0=2j3,cj+2,% =0 
as well as gq=a/3v. 

Accordingly, Eqs. (H.24), (H.25), (J.18) and (J.19) must be written as follows: 


bi=(1—@) —(@/3)Ao,0.0 , (K.1) 
b2=(1/2)(1—w) — (o/3)(Bo,o,o— A? 9 9/2) , (K.2) 
Ao,o,o=—(v+2 >, 19-9; . 

2>0 


ER Vs RI tN (K.3) 
Dea dy has sak; t 


+2H RL TN ys, citys at Gh; )/v®. (K.4) 


j,k>0 


By using the new variables introduced at the end of Sec. J, the summations 
in Eq. (K.3) for As,o,0 may obviously be given as follows: 


V4.0, = DG, (K.5) 
j>0 j>0 

DTS TU jan cet Da His caer « (K.6) 
2>0 3>0 


For carrying out the summations in Eq. (K.4) for Bo,o.o, the following 
consideration on the feature of Yang-Li-Hill lattice is very useful: On this 
lattice, once two separate chain polymer molecules approach each other to 
occupy a triangular cell of the lattice in common, the other parts of the two 
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polymer molecules can never approach each other to occupy another triangular 
cell in common. 


It follows, therefore, that the summations may be factorized as follows: 
3 hp; ete = (Da 2) (THOT gy.) 
j>0 k>0 


g,k>0 
ee) Sh) (K.7) 
j>0 k>0 


Pane are | 
NE Iss i 080k vo 
j,k>0 


=( う > I ; Bee) (os Ore ONY so) 
j>0 k>0 


=( 21 >.<+1 ( 2 Se) . (K.8) 
Since, 
§>,<G+1=Y i>, < シー 2>5>,<j41>,<j+20 
— £5j-1>,<j>,<jtic » 

it follows from Eq. (H.3) that 

25> .<3-1=[(@/2e) — (BP /3») Jo ’ 
where 

a=1=(1/p) , b=1—(/d), 
and 


B= (SBi(P—2) . (K.9 


Furthermore, since 
Ci =2j5—2(N5>,< 541 t95-1>.<5 
+ 255>,<j41>,<j42c $25j-1>,<j>,<s41c F25j-2>,<5-1>,<ic) » 
it follows from Eqs. (H.3) and (K.9) that 
21 55=[1— alo) + NS) lv é (K.10) 


When Eqs. (K.1) to (K.10) are combined, it follows finally that 


bi=1/p, 
b2=(1/2) + [1 —3(a/a) + (3/2)(a/)"] 
) —oH{1—2(a/a) +(1/3) (B/»)]- 
-[{1—2(a/) + (1/3) (B'6/»)} 
+2H{(1/2) (a/)—(1/3) (8P/»)}] . : (K.11) 
When athermal, i.e., H=1, this equation may be reduced to 


be=(1/2)—(1/2) (a?/a) 
—(1/3) (8b/v)[a—(1/3) (BPbo/v)] , (Athermal) (K.12) 
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and further to Eq. (1.6), by putting B=; =0. 


L. Second (Kramers-Wannier-Kikuchi) Approximation 


The second, i.e., Kramers-Wannier-Kikuchi, approximation method corre- 
sponds to the special case q=1 of the preceding formulation. In this 
approximation method, the calculation of osmotic virial coefficients is extremely 
long. Hence, the calculation of the first virial coefficient is given in the first 
place. Furthermore, since the calculation of the second virial coefficient is 
very much complicated, it is given only for the case of a two-dimensional 
triangular lattice in order to illustrate the method of calculation. Full ex- 
planation of the present method will be given in a separate paper including 
some applications to three-dimensional lattices. 


First Osmotic Virial Coefficient 


The first virial coefficient may be obtained by combining Eqs. (J.16) and 
(J.18) with Eq. (H.24). The summations appearing in Eq. (J.16) for ao. and 
Eq. (J.18) for Ao,o,o may be written by using new variables introduced at the 
end of Sec. J as follows: 


il Ons = DIGG, (LD 
j>0 j>0 

3 SATO GS eset sae Dy Md. < 944 » Cha) 
j>0 32>0 

SH ed yj. cita: t= | . (L.3) 
j>o0 j>0 


Since, in the present case, 


Ni>.<7t1=Yi>,<j4+1—2d5>, <j41>,<j+2c — 2 aj-1>,<j>,<j+ic > 


Nj>,cj+2=Yjr.cj+a— 2aj>,<jt+1>,<s+ac >» 
it follows from Eq. (H. 3) that 
2 i> .<it1=[(a/2m)—(B/3v)I0 , 
Baera= (BD/20)— (Bb/6Dv , 
where, as in Eq. (K.9), 


a=1—(1/p), b=1—(2/p), 


and 


B= (2185) (P—2) . (L.4) 


Furthermore, since 
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C5 = 05-207 §> << 54a bN5—1>,<5 +Ni9,c549tNj=a>,c4 


eS <j+1>, <ptac 253-15, <§>,<j4-1e + 2>9-a>,<jni>,<ic) ; 


it follows from Eqs. (H.3) and (L.4) that 
ZS =[L— Zale) — (28 b/w) + (BAND . (L.5) 


When Eqs. (L.1) to (L.5) are combined with Eqs. (J.16) and (J.18), it 
follows that 


do ,0= —[2—(a/o)—(B™b/a)] , (L.6) 
Ao,0,0o= —[3—(3a/@)—(38'b/w) +(BB/v)]Jv . «L.7) 


By introducing these equations into Eq. (H.24), the first virial coefficient 
may be found to be 


bi=1/p , (L.8) 


which may legitimate the foregoing calculation. 


Second Osmotic Virial Coefficient (Two-dimensional Lattice) 


As mentioned in the earlier part of this section, the practical calculation 
of the second virial coefficient is made for a two-dimensional triangular lattice, 
w=3 and v=2, for which Eq. (H.25) is simplified, after introducing Eqs. (J.17) 
and (J.19) into Eq. (H.25), as 


be=(1/2)—3H Sy gj ; ote ;el0®@—(1/2)(345,0—2Ad0.0) , (L.9) 
J 


leaving only the summation over clusters /, R, 0. 

In order to carry out this summation, the following consideration on the 
feature of Kramers-Wannier-Kikuchi lattice is very useful: On this lattice, 
there may be found sequences of triangular cells connected side by side on 
their lattice bonds. It may be seen, therefore, that if two separate chain 
polymer molecules approach each other to occupy a triangular cell of the lattice 
in common, there remains still another possibility that the other parts of the 
two molecules might approach each other to occupy another triangular cell in 
common. Hence, if this summation is written as ( 2a ban =, be) asin Eq. (K.7) 


of the previous section, there may be included impossible pairs of segments 

j;s and k;t of two separate polymer molecules of which two branch chains, 

starting j7;s and k;t, respectively, would overlap at a certain lattice point 

outside the triangular cell on which the pair j;s and k;¢ have been found. 
It follows, therefore, that the summation must be written as follows; 

Dd Fes. ee eC 2 bs) ae" ’ (L.10) 


j,k>0 
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where « is a correction factor for the aforementioned overlapping. « is called 
the overlapping correction factor, which may be expressed quite formally as 


follows: 
p 
kab ike for p-mer, 
Ke ate Ll) 
Tr 
ks is called the s-mer correction factor, since correction factors up to xs should 
be taken into account so as to correct the overlapping of s-mers. The mode 


of overlapping of two branch chains corresponding to «s'” is illustrated 


schematically in Table 3. 


Table 3. 


@+@ Indicates the segments are connected with either chemical or inner 
correlation bond. 


(a) Cluster Variables for Estimations of Overlapping Correction Factors 


For the estimations of kx» and xs, new variables given in Table 4 are in- 
troduced. 
With these new variables, «2 and ks may be given as follows: 


k= C5")*/0? , 
w= DOP LX (cgi ty cin 
+P esa types) , 
Ks” =2{> (NF part Fat at NHC} 42)}° 
LD pt to a+ 1} 49) }/v? . (L.12) 
The estimation of x:, s>3, may be performed in exactly similar manner 
by introducing cluster variables representing probabilities of finding larger 
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portion of polymer chain. Table 4. 


However, in the present paper, 
Van, (0,1) 


iy +t 


vy (1,0) 
Np. 


Vo Sow 


5) (+2 


A 


~ 5) +2 


the final results will be given 
for some simpler cases. 

For the evaluation of the 
first two summations in Eq. 


ae enw ais was se a oom fd 


(L.12), consider, on the original 
triangular lattice, a triangular 
cell, say abc, of which a lattice 
point, say a, has been occupied 
by asegment j. Furthermore, 
consider a group of six trian- 


(0 ,0)(1) (1)(0,0) 
Ny a+ VE 


(1,0)(0) (00,1) 
pith Ds. tl 


(0,0)1) (DO,0) 


( (1,0)(0) (0)(0, 1) 
25). G+2 75) .G+2 TG 


J). +2 2) G+2 


gular cells surrounding the @- Indicates the lattice bond is occupied by either 
lattice point a, of which one is orange 
the triangular cell abc. On this group of six triangular cells, there will be found 
various modes of configurations, with respect to the triangular cell abc, of two 
branch chains j, 7—1, j—2,--- and j,j+1,7+2,---, starting the segment 7. Each 
of the modes of configurations is specified by index a. When the configuration 
of the two branch chains on this group of six triangular cells is found in the 
a-th mode, the segment 7 is written as 7; a. The subindexing a corresponds to 
a coarse-grained grouping of the subindexing s used hitherto. The probability of 
finding a segment j:a@ on a lattice point of a triangular cell is written as 4; , a. 

By the restriction imposed on the occurrence of intramolecular correlations 
at the beginning of the present formulation, the subindexing a may be made 
by observing whether each of the six lattice bonds, ab, ac, etc., diverging from 
the lattice point a, is occupied or not by either chemical or inner correlation 
bond. 

The probability of finding a certain mode of occupancy of the six lattice 
bonds by chemical and inner correlation bonds is determined by shape factors 
Bn, BY, Bl, and BY. Hence, x;,« for the present formulation may be given 
in terms of v and these four shape factors after elementary but laborious 
calculation. 

By summing %;.a over all of a which are compatible with €; and €;°.”, 
respectively, and also by summing them over j>0, the first two summations 
in the equations for «2 and xs of Eq. (L.12) may be obtained. The final 
results will be given shortly. 

Quite similarly, by considering the configuration of two branch chains, one 
starting the segment j and the other starting the segment j+1, (+2), of 
cluster j>,<j+1, (j>,Cj+2), on a group of ten triangular cells surrounding 
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the cluster j>,<j+1, (j>, Cj+2), the summations over 7 in Eq. (L.12) may 
be obtained. 


(b) Second Coefficient for Tetramer 


For tetramer, p=4, the summations appearing in Eq. (L.12) may be given 
as follows: 

DiGi = (1/24) {Ba + fa" + (1— Ba) +(1— Ba) +282 85 
+ Bo(1— Bs) + Bs(1— Bo") + Bo") Bs“(1 — Bo) 

+ a! Ba'(1—Ba®)-+(1— Ba) o((1— Bo) 
+(1—fs'?) Be" (1 — Bs") + B2'(1— Bo")(1 — Bs) 
+ Pa'(1—2)(1— Bo) 
+2(1—Be')(1—fa')(1— Ba) 

+2(1—Ps'®)(1— Be')(1—Ba')}o , 

SG: 1) =(1/24){ Bo" + Bs") +(1— 2") +(1— fs) +282") 83 
+ Be (1 — Bs) + Bs(1— Bo") 

+ Ba'*(1—Bo')(1— Bo!) + Bo — a)(1—Aa')}o , 

7 eM 5 4g =(1/48){2 Bs! Bo (1 — Bs‘) + Ba (1— BD) (1— Bs") 
+(1— Bs") B21 — Bs") 

+ Bo" Bs\(1— BD (1— Bs"))(1 — Bs‘) 
+ P2'(1—Pa')(1—Ba)\(1—a)}o , 

Dj = (1/48){ 280 pa D1 — Bol) + Ba C1 — Bo )(1—Bs™) 
+(1—Bs'®)Bs""(1— fa) 

+ Be") Bs'2)(1— Bo"))(1 — Bs“)(1— Bs") 
+ Bo!(1— Ba ®)(1—Ba'9\(1—o)}0 , 

DRL; ra = (Pa? /48{(1— Ba). —Pa®)}0 , 

Da = (aD/48 1 — Ba! \(1—Ba')}o 

BY 5 = (1/48) {Ba Bo — Ba") +1— Ba a1 — Ba") 
+(1—fal®)Ba'(1— Ba) 
+2(1—Ba!®)(1—Ba')(1— fs) 

Ba (1— Bo \(1—Bs'X1.— Ba") 
+2(1—Pa')(1—a®)(1—Bs')(1—Bo)}o , 

DR = (48 Bo Be' (1 — Bal) +-(1 — Ba!) BoC — Ba) 
+(1— Ba) Bs%(1 — a) 
+2(1—Bs')(1—Ba'?)(1— Ba") 

+ Ba®(1—Ba®)(1—Ba)(1— fa") 
+2(1—fa!®)(1—Ba)(1— Bs" \(1—B)}0 , 
DN, = (Bs'"/48 1 —Bo™)} 0, 
ra = (N48 {B41 — Ba , 
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DI 1 =U/48) (Bs aD 1 — oD) + Bal D1 — Ba) 
+(1— BaD) a1— Bo") + Ba'(1— oN pa) 
+ Bo!” Bs'(1 — BaD 1 BY , 
BR = (1/48) { Be!” BaD — a) + fC — Ba) 
+(1— Bu!) Bs(1— fo) + Ba'(1 aD Batt) 
+8 Ba'™(1—Bo®)(1— Ba®)(1—Bs"™)}o , 
DI sg =(Pa? /48){(1— Ba). — Ba} , 
D790), 4 = (62/481 Ba)\(1— Pa) . (L.13) 


With these equations, the overlapping correction factors kz and xs may be 
calculated explicitly by Eq. (L.12). «s is evaluated in a similar manner in the 
case of a rigid polymer chain and is given in Table 5. In the case of a 
flexible polymer chain, « is ignored. With these correction factors, the second 
virial coefficient for tetramer may be calculated by combining Eq. (L.9) with 
Eqs. (L.6), (L.7), (L.10), (L.5) and (L.11). The numerical results for athermal 
solution, H=1, are summarized in Table 5. 

The rigorous values of the second virial coefficient are given in the last 
row of the table which are obtained by computing geometrically the excluded 
volume of polymer chain. The details of the method have been given in Sec. 
H of Ref. 13). 

For rigid tetramer chain, there may be found perfect accordance between 
the calculated and rigorous values of the second virial coefficient. However, 
in the case of flexible tetramer chain, since the tetramer correction factor has 
been ignored, the calculation gives an under-estimated value. Nevertheless, 
the calculated value seems to be correct within an error of 2%. 


(c) Second Coefficient for Extremely Long Chain Polymer 


In the case of extremely long chain polymer, fp-oo, rigid shape chains 
are of little interest. Hence, shape factors are chosen in advance as 


Bi =Bay ’ Bj = Be) ’ 
and 
= ; tou all 7. (L.14) 
The summations appearing in Eq. (L.12) may be given, after making 
p-o, as follows: 


DiGi =(1/6){ 64, +284,(1— Bay) +381 — Bin)? 
+(2—Br))(1—Bu)*}o , 
dCi) = (1/6) (262, (1— Bay) + Bay 1 — Bay)? 
+ Bra(1—Bay)?}o , 
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Table 5 
Rigid* 
。 () = R31) =0, (=f) =0, Bo) =f, =0, | Bot) =BsD=0, 
Ba) = pitt) =O, Ra 0 の ご 1, B® = Ba?) =1, Bo") = B32) =1, 
B33) =0 Bs(3) =0 B3(3)=0 Bs =1 
(a) (b) (ec) (d) 
(Cv (12/24)? (12/24)? (12/24)? (12/24)? 
Ke (6/24)? (5/24)? (4/24)? (4/24)? 
- 0 ks(11/96 sD1/36 ra(D1/72 
3 
ka(2)1/384 
ka 2)1/1152 ka D1/144 
Ka 0 ks)3/1152 0 ka 1/288 
Fa)1/288 
See ie 3/16 73/384 7/36 7/36 
—,0 7/4 7/4 7/4 7/4 
— Ao,o.0 9/4 9/4 9/4 9/4 
be 
cale. 
156/384 153/384 148/384 148/384 
be 
rigorous 


* es, (b2 rigorous 116/384) is omitted, since it violates the restriction on intra- 
**  (a)~(g) occur in the following mole fractions: (a) 2/18, (6) 4/18, (c) 1/18, (d) 1/18, 
*** Weighted mean of be for (a)~(g). 


DY) aera PSL) ar 
=(1/24)({26%, (1—Bin) + (3 +28c2)(1—Bin)*B wy 
+2(1+ Brey) By(1— Bi 
+28(2)(1—Bi2))(1—Bay)*}o , 
Pa) Laer a> >.) ree 
= (Bi1)/12){ Bt, (1—Bay) +(2— Bi) 62, (1— Bay)? 
+ Biy(1— Bay)? + Bball— Be) (1— By) Fo , 
Binh) , = TV 
= (1/24) {(8—2B.2)) Bt, (1— Bay)? 
+(6—5Br2))Bi(1— Bay)? 
+2(2—3B (2) + Bt )(1— Bin) Fv , 


sc 
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Tetramer (p=4) 


Solutions 81 


Flexible 
pt 21, '9,07=0, Bo) =1, 8 w=0, gp wag Wt, Rie Wa] 
Bo?) = Bg?) —0, Bo") =0, 22 =1. Ba) = B®) =0, ie te 
Bs =0 83°) =0 B33) =0 Ba =1/2 | 

@ (f) (9) a 
he / Perfectly 
ee dvs flexible** 
(11/2492 (11/2492 (10/24)2 (17/36): 
(4/24)2 (4/24)? (4/24)? (5/27)2 
ka()1/96 a(t) /96 4 a1420/(36. 26) 
ks (21/288 ks 2)1/288 ks (D143/(36.26) 
0 0 0 agnored 
97/576 97/576 7/48 2747 (35-28) 
5/3 5/3 19/12 61/36 
17/8 17/8 2 78/36 
~137.4/384 
132/384 132/384 116/384 
~139.3/384** 


molecular correlations. 


(e) 4/18, (f) 4/18, (g) 2/18. 


(0,0) (1) — 20 (0,0) 
Dy aay ae aa UE Cea Pe 


=(Bi1/12){ Bt, + B3,,(1— Bay) +282) 8%, (1 — Bay)? 


+ Bia) Biy(1— Bay)? + 4.1 — Bay) to , 

bay why rei eH publ are 
= (1/24)426, (1— Bu) +21 + Bee) 62, 1— Ban? 
+382) Ba (1— Bay) + 83 (1— Buy)! bo , 


(1,0)(0) — (0) (0,1) 
EA = QU peg +8 


= (Bi/12){ 6%, (1— Buy)? +(1— Bra») Bay — Bay)*}o « 


(L.15) 


With these equations, x2 and kxs may be calculated by Eq. (L.12). ts, s>3, 
are all ignored. The numerical results for athermal solution, H=1, are 
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Table 6 Extremely long chain polymer (p>) 
Rigid : Flexible 
=0,8:020,5 li pimtle 0 bi o> Aree mae = 
psi Os, " ear r B(3)=1/2 Bis) =1/2 
i) i perfectly flexible 
paar pO VAS fa 
(ZC (6/18)? (3/18)2 (6/18)? (5/18)? 
Ke (6/18)? (3/18)? (3/18)? (29/162)2 
k 0 0 ks( 1/36 ga)232/(3.94) 
: x3(2)980/(3-95) 
ks(s>3) 0 0 ignored ignored 
See 0 0 1/18 4924/(3-98) 
—Q0,0 5/3 4/3 5/3 14/9 
— Ao,o,0 2 3/2 2 11/6 
be 2/12 34979/ oe 95) 
cale 4/12 1/12 =0.167 =0. 
_ be 0.333 0.083 be 
rigorous ee oS” os. 
1/4 
wee «=| = 0.250 
3/16 
Sood a =0.188 
1/9 
0.111 


summarized in Table 6. 
As seen from Eq. (L.14), two cases, Bi) =f2)=0 and Bu) =1, Ba, =0, correspond 
to a rigid stretched chain and a zigzag chain, respectively. The results for 


these two special cases are included in the table. 


The case Bu,=1/3, Bs,=1/2 


and Bs,=1/2 corresponds to a perfectly flexible chain. Since x, for s>3 has 
been ignored, the value of the second virial coefficient for a flexible chain is 
under-estimated. However, the calculated value seems to be quite reasonable. 
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M. Feature of Higher Approximation 


In the preceding sections, the second osmotic virial coefficients bz for the 
first, (1+@/3v)-th and second approximation methods have been calculated as 
three special cases of the (1+gq)-th approximation. In the first approximation 
method, it has been impossible to take the shape of a polymer chain into 
consideration. In the (1+/3yv)-th approximation method, the shape has been 
taken into consideration only partially through the shape factor Brthers 
more, in the second approximation method given above, the shape has been 
expressed by three kinds of shape factors 68;0, 8B;% and B;3, in terms of which 
rigorous values of bz have been calculated for all but one of eight kinds of 
possible shapes of a tetramer as shown in Table 5. However, when chain 
polymer molecules are of rigid linear shape and accordingly §;‘”=£;'=£;'*)=0, 
it may easily be shown that both of the expressions of bz for the (1+/3v)-th 
and second approximation methods reduce to the expression for the first ap- 
proximation method. In other words, the shape of a polymer chain utterly 
loses its meaning in the first approximation method. 

Calculations may give, for a two-dimensional triangular lattice, the follow- 
ing values of b2"'+7] for an extremely long chain polymer (p=%), as Bj;W 


increases from zero to unity, 


lst approximation: bLN=4/12 , 
(1+0/3v)-th approximation: b{1+ ©/8vJ—4/12-33/12 , 
2nd approximation: b!2I=4/12—1/12 , 


where it must be noted that b2=4/12 and 1/12 are the respective rigorous values 
as shown in Table 6. Thus the second approximation method provides a rigorous 
method for the estimation of the second virial coefficient at least in the two 
limiting cases of §;‘=0 and 1 for all 7. This is because of the fact that, in 
the two cases, there can never be found any intramolecular correlation between 
two segments j and j+s with s=3 of which the occurrence has been prohibited 
by the restriction premised at the beginning of the present formulation. 

It may be suitable here to return to the consideration of the nature of 
this restriction in detail. Suppose there is a chain polymer molecule on a two- 
dimensional triangular lattice. The number of different kinds of shapes, which 
may be taken by the polymer chain, is three for a trimer (p=3), eight for a 
tetramer (p=4), thirty two for a pentamer ( p=5) and about one hundred and 
twenty for a hexamer ( p=6) of which none, one, eight and forty six, respec- 
tively, violate the restriction. In brief, the defect of the present formulation 
becomes more marked as the length of a polymer chain, i.e., p, increases. 
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In order to improve further the present theory, the moderation of this 
restriction may be useful. For this, natures of various kinds of intramolecular 
correlations must be considered. When an intramolecular correlation between 
two segments j and j+s of a single polymer molecule is found on a lattice 
bond, say ab, of a two-dimensional triangular lattice, 7 on a and j+s on 3B, 
there may be found that the portion of the polymer chain between the segments 
j and j+s, together with the lattice bond ab, encloses a certain region of the 
triangular lattice. Next, consider a sheath enveloping this closed region. If 
there is found none of hexagonal or polygonal sequences of lattice points within 
this sheath, the intramolecular correlation between the segments j and j+s 
is called an intramolecular correlation with the first type. On the contrary, if 
any sequence is found, it is called an intramolecular correlation with the second 
type. 

The intramolecular correlation with the first type may be taken into account 
within the framework of the second approximation method by merely increasing 
the number of cluster variables by the number of higher intramolecular cor- 
relations added. On the contrary, the intramolecular correlations with the 
second type can never be taken into consideration by the second approximation 
method, since such a configuration of the chain on hexagonal or polygonal 
sequence of lattice points cannot be transcribed from the original triangular 
lattice onto the corresponding Kramers-Wannier-Kikuchi lattice. 

It may be seen, therefore, that the prohibition of the intramolecular cor- 
relations with the first type is made only for the convenience of facilitating 
the formulation of the second approximation method, while the neglection of 
the intramolecular correlations with the second type is intrinsic in this method. 

The removal of the restriction on the occurrences of the intramolecular 
correlations with the first type makes it possible to decrease the number of 
shapes which have been omitted from the present formulation from one to zero 
in the case of a tetramer, from eight to zero in the case of a pentamer and 
from forty six to only one in the case of a hexamer. However, in order to 
improve further the present theory by removing the restriction on the occur- 
rences of the intramolecular correlations with the second type, it is necessary 
to proceed to higher approximation methods, in which hexagonal or polygonal 
correlations of lattice points should be taken into account. 
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IV. Summary 


In the present paper, a systematic method of constructing combinatory 
factors has been developed, and its application to chain polymer solution has 
been made. So long as the second osmotic virial coefficients of globular polymer 
solutions are concerned, the gas-like model theory has provided an elegant 
method for rigid polymers having geometrical shapes expressible analytically. 
However, the situation is somewhat different in the case of chain polymer 
solutions where the gas-like model theory has been confronted with mathema- 
tical difficulties. 

The lattice model theory, to which some contribution has been made in 
the present paper, is effective in giving concrete results when applied to chain 
polymer solutions. The present author has shown in this paper that it is 
possible to take into account simultaneously the size and shape of a polymer 
chain through p and f;"’, and the temperature and the heat of solution through 
H by higher approximation methods. 

The lattice model theory has involved some ambiguity in the physical 
nature of combinatory factor. Such ambiguity has been removed by introduc- 
ing “Bethe”, “ Yang-Li-Hill ” and “ Kramers-Wannier-Kikuchi” lattices. 

The lattice model theory has another advantage of giving expressions over 
the whole range of concentrations, while the gas-like model theory makes its 
approaches only from higher dilutions. However, the present author has 
confined himself to the theoretical calculation of the second osmotic virial 
coefficient in order to investigate the feature of the lattice model theory clearly. 

In the improved lattice model theory given in Sec. L, the second osmotic 
virial coefficient is expressed analytically in terms of shape factors, and calcula- 
tions give rigorous values of the second osmotic virial coefficient for a two- 
dimensional triangular lattice in a large number of cases. The same method 
may be applied to other types of lattice. However, in order to explain the 
principle of this general method, the formulation has beer limited to a two- 
dimensional triangular lattice. Further applications, including those to three 
dimensional lattices, will be given in a separate paper. 


(Contributed Feb. 28, 1961) 
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